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ABSTRACT 


The  research  efforts  for  the  first  year  on  Project  Themis  at 
Texas  A&M  University  have  been  di  'acted  along  five  primary  research 
paths:  analysis  and  prediction  of  atmospheric  st:  icture  through  use 
of  dual-frequency  radars  and  associated  ground  facilities;  cumulus 
cloud  dynamics;  modeling  studies  of  the  frictional  boundary  layer; 
Installation  of  two  automatic  meteorological  collection  stations  in 
an  Instrumented  watershed  which  is  in  the  pattern  of  the  radars;  and 
design  of  a  facility  for  direct  measurement  of  the  fluxes  of  momentum, 
heat,  and  vapor.  The  activities  under  these  various  headings  are 
reviewed  separately  herein. 
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INTRODUCTION 


The  meteorological  research  program  under  Project  Themis  at 
Texas  A&M  has,  since  its  inception,  been  divided  into  four  main 
sub-areas  as  listed  below. 

Sub-task  a.  Cumulus  Cloud  Development  and  Dissipation 

Sub-task  b.  Analysis  and  Prediction  of  Atmospheric  Structure 
Through  Dual  Frequency  Radar 

Sub-task  c.  Analog  Modeling 

Sub-task  d.  Turbulent  Flow  Definition  Through  Reynolds 

Stresses  Measurement 

The  single  research  goal  under  this  program  is  the  development 
of  meteorological  -  mathematical  models  to  improve  objective  forecast 
capability  on  all  scales.  Of  course,  it  is  an  oversimplification  to 
group  research  on  all  atmospheric  scales  under  a  single  heading  and 
therein  really  lies  the  reason  for  the  subdivision.  However,  some 
need  exists  to  show  how  these  sub-areas  make  up  the  whole  and  it 
is  to  this  point  that  this  introductory  section  is  directed. 

The  Themis  research  program  derives  in  part  from  earlier 
research,  conducted  for  the  Atmospheric  Sciences  Laboratory  of  the 
United  States  Army,  which  led  to  the  development  of  a  model  of 
the  atmospheric  boundary  layer.  The  specific  application  of 
this  model  is  from  the  surface  to  an  elevation  of  1000  meters 
and  from  the  surface  to  a  depth  of  two  meters.  This  model 
and  associated  research  revealed  several  deficiencies  in 
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fundamental  meteorological  understanditg  that  comprise  a  pert  of 
the  investigations  pursued  in  the  sub-arees  under  Project  Themis. 

The  analog  modeling  portion  (Sub-task  c)  is  not  a  continuation 
of  this  earlier  modeling  study  but  rather  an  investigation  by  analog 
means  of  the  fundamental  gaps  in  meteorological  understanding, 
par  '.cularly,  on  the  meso-scale.  A  significant  example  in  this  vein 
is  the  diabatic  profile  of  wind,  temperature,  or  vapor  pressure. 

Of  course,  several  "profile  laws"  exist  and  a  couple  of  these  have 
been  cited  so  often  in  the  literature  that  they  enjoy  the  status 
of  being  widely  known  and,  unfortunately,  widely  accepted.  Un¬ 
fortunately,  since  they  are  impractical  in  the  sense  of  theoretical 
derivation  and  physical  realism  except  for  rather  narrow  ard  limited 
applications.  Another  example  concerns  the  exchange  coefficients, 
which,  though  widely  used,  are,  as  yet,  without  functional  definition 
and  understanding  is  at  best  confused.  Historically,  as  well  as 
because  of  some  practical  necessities,  this  definition  has  been 
sought  through  mean-flow  characteristics  since  the  latter  are 
the  only  parameters  that  can  be  routinely  measured.  To  date, 
equivalence  has  not  been  determined  and  simultaneous  direct 
measurement  of  Reynolds  Stresses  and  pertinent  mean  parameters, 
under  Sub-task  d,  is  a  logical  path  to  follow. 

A  deficiency  of  the  noted  model,  and  other  models,  rests 
in  the  failure  to  include  cloud  formation  and  dissipation. 

For  limited  model  areas,  this  is  not  too  serious  a  problem  but, 
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It  be com s  a  aoat  serious  problem  aa  the  forecast  area  Is  extended 
vertically  and  horisontally,  and  the  cumulus  dynamics  studies 
pursued  under  Sub-task  a  are  addressed  in  part  to  this  problem. 

Of  course,  turbulent  flow  definition  under  Sub-task  d  is  also 
important  to  the  cumulus  dynamics  study. 

In  general,  existing  models  do  not  Include  precipitation, 
except  in  the  total  water  balance  sense,  nor  incorporate  passage  of 
frontal  systems.  The  reason  for  this  lack  is  twofold:  1)  the 
functional  definitions  of  involved  phenomena  are  unknown;  and  2) 
the  measurement  capability  to  provide  necessary  model  input, 
assuming  such  understanding  existed,  is  not  routinely  available. 

This  introduces  the  dual  frequency  radar  studies  under  Sub-task  b 
which  seeks  to  provide  possible  input  through  radar  surveillance 
as  well  as  permit  Incorporation  of  precipitation  modeling.  These 
studies  also  supplement  the  cumulus  dynamics  modeling. 

As  noted,  these  brief  introductory  remarks  are  intended  only  to 
show  the  relationship  of  the  various  sub-tasks  under  Project  Themis 
to  the  overall  project  goals.  Detailed  information  under  these 
activities  are  reported  separately  on  the  following  pages. 


II.  CUMULUS  CLOUD  DEVELOPMENT  AND  DISSIPATION 
K.  C.  Brundidge 


A.  Introduction 

The  long-range  objectives  of  this  portion  of  the  Themis  program 
(Sub* task  a)  are: 

1)  to  conduct  numerical  experiments  on  cumulus  cloud  initiation 
In  which  various  conditions  may  be  asalgned  with  respect  to 
heat  and  moisture  sources,  ambient  wind,  and  exchange  coef¬ 
ficients, 

2)  to  extend,  where  possible,  previous  theoretical  models  of 
cumulus  cloud  development, 

3)  to  compare  the  computed  cumulus  development  with  available 
observational  data. 

The  accomplishment  of  these  objectives  will  depend  in  part  upon  the  ob¬ 
servational  program  conducted  in  Sub-task  b,  and  the  determination  of 
exchange  coefficients  as  carried  out  in  Sub-task  d. 

The  specific  goal  of  Sub-task  a  for  this  first  year  has  been  to 
formulate  the  mathematical  procedures  and  models  to  be  used  in  the  study. 
The  preparation  and  testing  of  the  initial  computer  programs  for  the 
various  models  was  then  to  be  accomplished  during  the  next  18-month 
period.  In  actuality,  there  has  been  some  blending  and  realignment  of 
these  two  phases  during  this  first  year.  As  a  result,  the  moist  con¬ 
vection  model  has  not  been  completely  formulated  as  yet;  however,  an 
initial  program  for  a  two-dimensional  dry  convection  model  has  been 
written  and  has  undergone  some  testing.  This  is  in  keeping  with  our 
decision  to  concentrate  on  the  simplest  of  the  various  possible  models 
as  a  first  step.  In  particular,  we  have  devoted  most  of  our  attention 
to  a  model  of  two-dimensional  dry  convection  with  zero  ambient  wind 


field.  The  advantages  of  this  procedure  are  that  we  cpxf check  the 
accuracy  of  our  computer  program  with  considerably  less  computer  time 
than  would  be  required  for  a  more  complicated  model,  and  many  of  the 
developed  sub-routines  will  have  applicability  in  a  more  complicated 
model.  Furthermore,  we  can  determine  the  influences  of  various  initial 
conditions  and  heating  rates,  multiple  heat  sources,  and  various  assump¬ 
tions  for  the  eddy  exchange  coefficients  on  convection  before  bringing 
in  the  added  influences  of  latent  heat  release. 

In  the  following  sections,  we  shall  present  the  experimentation 
program  which  has  been  set  up  as  a  guideline,  the  details  of  the  various 
models  which  are  to  be  studied,  the  problems  that  have  been  encountered 
and  their  solution,  and  the  numerical  procedures  which  have  been  inves¬ 
tigated. 

The  basic  equations.  The  basic  set  of  equations  being  used  in  this 
study  is  the  anelastic  system  developed  by  Ogura  and  Phillips  (1962). 
This  system  permits  gravity  wave  solutions  but  filters  out  compression 
waves.  The  formulation  of  this  system  for  moist  convection  will  be 


stated 

at  a  later  point. 

For  dry  convection  and  in  dimensional 

form. 

it  is 

iH  _  _r 
dt 

®-|x  +  KhJ^V®  u  +  a(z)uz  J 

1) 

—  -  -C. 
dt  ^ 

®|~  +  K*|V  v  +  a(z)vz  ] 

2) 

dt  “ 

®li  +  8I  +  W  +  a(r)w?] 

3) 

5 


v-pof  ■  o 

e'  +  -(2)0!'] 

where  p0  "  ^  (l  ■  p*) 


4) 

5) 


X  -  1 

i  K 


and  a(z)  -  Po"1  ,  0  - 

P  is  a  reference  pressure  (say,  1,000  mb),  8  it  a  constant  reference 
potential  temperature,  K„  and  are  the  eddy  coefficients  for  momentum 

it 

and  heat  exchange,  respectively,  ^  ™  p  where  p  is  pressure,  and 
tt'  is  the  deviation  of  tt  from  the  distribution  tt0  =  1  -  gz,  i.e., 

Tf'  "  TT  -  TT0 

6'  is  the  deviation  of  potential  temperature  from  the  value  0,  l.e., 

e'  -  e  -  e 


and  other  symbols  in  1)  -  5)  have  the  usual  meteorological  meaning. 

This  system  has  been  used  by  Ogura  (1963)  and  Orville  (1964,  1965) , 
among  others,  in  studies  of  convection  in  a  vertical  plane  or  in  an 
axially-aymmetric  cylindrical  coordinate  system.  These  studies  have 
been  restricted  to  shallow  convection,  i.e.,  where  the  top  of  the  model 
is  at  about  3  km.  For  the  latter  situation,  Ogura  and  Phillips  (op.  cit.) 
have  shown  that  4)  may  be  replaced  with  V°C  ■  0.  Then,  the  equations 
amount  to  the  incorpressible  Boussinesq  approximation.  With  this 
combination  of  circumstances,  one  may  make  use  of  the  streamfunction  and, 
furthermore,  may  replace  1)  and  3),  or  their  cylindrical  c>  uiiiates 
equivalents,  with  the  vorticity  equation.  Consequently,  it  is  never 
necessary  to  compute  tt'. 
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However,  it  is  our  desire  eventually  to  consider  three-dimensional 
regions  and  deep  convection  (top  of  the  region  at  about  10  km.).  for 
this  purpose,  it  is  necessary  to  work  with  the  primitive  equations  as 
given  by  1)  -  5)  and  this  requires  a  means  whereby  tt'  can  be  deter¬ 
mined.  This  is  especially  critical  for  moist,  deep  convection  for  then 
n'  must  be  known  in  order  to  compute  the  saturation  mixing  ratio.  A 
diagnostic  equation  for  computing  tt/  is  obtained  by  forming  the  diver¬ 
gence  equation  from  1)  -  3)  and  may  be  written  as 

Cp  0 \7  0  po  V  tt ,  «*  0^^“  +  V.poIP  -  '?®[Po$*^*l  6) 

AAA  AAA 

where  [F  *  K*  (1^0+  j^v+  kv®w)  +  a(z)  ^  (iu2  +  Jvz  +  kw2 ) . 
Equations  1)  -  6)  form  the  basic  system  with  which  we  are  working, 
aside  from  modifications  that  must  be  brought  in  when  an  ambient  wind 
field  is  considered  or  when  moist  convection  is  occuring. 

The  basic  model.  Natural  convection  in  the  atmosphere  is  a  mechanism 
whereby  heat  and  moisture  are  transported  upward  from  their  source 
region  a..  U.e  earth's  surface.  For  air  mass  type  convection,  the  lapse 
rate  may  be  superadiabatic  for  several  hundred  feet  above  the  ground. 
Accordingly,  we  have  chosen  to  work  with  models  which  incorporate  these 
features. 

We  consider  a  grid  system  in  which  Ax,  Ay,  and  Az  are  100  m, 
except  for  a  layer  resting  on  the  ground  for  which  Az  is  only  10  m. 

The  surface  temperature  is  a  specified  function  of  time  which  differs 
from  one  place  to  another  so  that  "hot  spots"  may  develop.  Heat  (and 
moisture)  is  transferred  from  the  ground  through  the  10-m  layer  by 
turbulent  diffusion. 
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The  initial  potential  temperature  distribution  which  has  been  chosen 
corresponds  to  conditions  in  the  Dallas ,  Texas  area  for  June  15,  1962  at 
0900  CST.  It  is  synthesised  from  the  radiosonde  observations  at  Anon 
Carter  Field,  observations  on  the  Cedar  Hill  tower  (Gerhardt  et  al,  1962), 
and  observations  on  a  mlcrometeorologlcal  tower  (Clayton,  1963)  located 
about  150  yds  from  the  Cedar  Hill  tower.  Since  we  have  not  as  yet 
entered  into  the  study  of  moist  convection,  the  moisture  distribution  has 
not  been  abstracted;  however,  the  data  are  at  hand  when  the  time  comes. 
Incidentally,  for  the  day  in  question  Love  Field,  Dallas  reported  cumulus 
clouds  with  bases  at  3000  ft  starting  at  1034  CST.  This  potential  tenq>er- 
ature  distribution  is  given  in  Table  II-l.  It  may  be  noted  that  the 
lapse  rate  is  superadiabatic  through  the  first  110  m. 


Table  II-l.  Initial  Temperature  Distribution 


The  experimental  program.  The  experimental  program  which  has  been  set 
up  is  given  below.  This  listing  does  not  necessarily  imply  the  order 
in  which  the  various  experiments  are  to  be  handled. 

Dry  Convection 

1.  Single  heat  source 

a.  Numerical  procedure  experimentation 

1.  Upstream  differencing 

2 ,  Others 

b.  Other  things  constant,  expand  X  dimension 

c.  Other  things  constant,  change  power  of  heat  source 

d.  Other  things  constant,  change  Initial  lapse  rate 

e.  Other  things  constant,  introduce  non-linear  diffusion 

f.  Other  things  constant,  introduce  ambient  wind  field 

1.  Vertical  shear 

2.  Vertical  and  horizontal  shear  (three-dimensional) 

g.  Combinations  as  seen  fit,  particularly  of  e.  and  f. 
ii.  Two  heat  sources 

a.  Other  things  constant,  vary  distance  between  sources 

b.  Other  things  constant,  vary  relative  power  of  sources 

c.  Other  things  constant,  vary  relative  size  of  sources 

d.  Other  things  constant,  introduce  non-linear  diffusion 

e.  Other  things  constant,  introduce  ambient  wind  field 

1.  Vertical  shear 

2.  Vertical  and  horizontal  shear  (three-dimensional) 

f.  Combinations  of  the  above  as  seen  appropriate 
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111.  Cumulus  Initiation 


a.  Single  heat  source,  surface  Inactive  as  moisture  source 
1..  All  Items  In  1 

2.  Other  things  constant,  vary  moisture  distribution  In  x  and  a 

b.  Single  hefct  source,  moisture- active  surface 

1.  All  items  in  ill  a 

2.  Other  things  constant,  vary  moisture  sources  at  ground 

c.  Two  heat  sources 

1.  Surface  Inactive  as  moisture  source 
(a)  All  items  In  111  b 
iv.  Cumulus  dissipation 

B.  The  Two-Dimensional  Dry  Convection  Model 

General  Features.  As  stated  previously,  this  model  has  received 
moat  of  our  attention  to  date.  It  is  the  most  simple  model  with  which 
to  begin,  although  it  Is  somewhat  unrealistic  In  that  It  assumes  that 
motion  occurs  only  In  the  x-z  plane. 

We  visualize  the  region  shown  in  Fig.  II-l.  It  is  assumed  thst 


Fig,  II-l  Two-Dimensional  Dry  Convection  Model  Area 
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initially  In  this  ragion  the  pressure  and  potential  temperature  surfaces 
are  horisontal  and  no  motion  exists.  At  the  ground,  differential  rates 
of  heating  take  place  so  that  In  a  short  time  the  stippled  area  In  Fig. 
11-1  becomes  a  hot  spot  relative  to  the  remainder  of  the  lower  boundary. 
Heat  diffuses  into  the  region  from  the  lower  boundary.  The  differential 
heating  then  acts  as  the  energy  source  to  initiate  and  maintain  the  mo¬ 
tion  field.  The  flow  will  be  symmetric  with  respect  to  the  vertical 
line  x  ■  0,  hence  we  need  only  consider  the  region  with  solid  boundaries 
in  Fig.  II- 1.  We  also  assume  symmetry  exists  across  the  vertical  x  *  L. 

For  testing  purposes,  we  are  letting  L  ■  3  km,  H  **  4  km,  and  the 
width  of  the  hot  spot  be  1  km.  The  temperature- change  rates  imposed  at 
the  bottom  boundary  are  as  shown  in  Table  II- 2. 

Table  II-2.  Surface  Temperature  Rates 
Time  [CST] 


09-10 

10-12 

12-14 

14-15 

lot  Spot 

2  deg  hr"1 

1  deg  hr"1 

0.5  deg  hr"1 

0.2  deg  hr"1 

Exterior  to 
lot  Spot 

0.8  deg  hr"1 

0.4  deg  hr"1 

0.2  deg  hr"1 

0.2  deg  hr"1 

System  of  equations.  The  system  of  equations  for  this  model  are  obtained 
directly  from  1)  -6).  They  are  as  follows: 

uT  -  -u  u*  -  wuz  -  Cp©n/  +  Km£uxx  +  u^  +  a(*)Uj  j  7) 

w,  -  -uwx  -  wwj  -  4  ®n/  +  ^  +  K*  [wxx  +  w„  +  a(z)wz  J  8) 
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] 


J 

| 

i 


i 


r 

i 

j 


po  (ux  +  Wj J  +  w(f.0)j  -  0  *) 

0t  -  -ue; •  we;  +  Knfe/x  +  e/,  +  a  <»)«,']  w 

whore  the  subscrlpCs  indicate  differentation.  Wien  the  continuity 
equation  9)  is  used  to  replace  teraa  on  the  right-hand  aide  of  the 
divergence  equation,  the  latter  take*  the  form 

Cj'ru  +  +  •T'i  ■  Cl(j9/  +  al'^)  Qtj-  <—•  ♦  ^  w,  -  CV  J 

+  Cll^U,*,  -  H 


where 


.  -  £  <p„>, 


01 

03  *^s 

_  8.75  C3l\ 


Ue  have  chosen  as  values  for  the  edl/  coefficients  i n  7),  1)  and 
10)  the  following:  «  K,  ■  7.5x10s  tW3  ^oc"1  at  Mic  10-e  level  but 

K„  -  K*  -  4xL0 4  col9  sec"1  at  all  higher  levels.  It.  »ay  be  noted  chat 
while  ^  and  are  not  independent  of  height,  this  fact  her*  been  ignored 
in  the  formulation  of  the  eddy  diffusion  term  in  V; ,  8)  and  10) . 
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Boundary  conditions.  The  boundary  conditions  were  chosen  to  provide  the 


symmetry  conditions  mentioned  previously  at  the  lateral  boundaries, 
zero  motion  at  the  upper  and  lower  boundaries,  slip  flow  at  the  top  and 
side  boundaries,  and  a  thermally' active  lower  boundary.  These  conditions 
are  as  follows: 


Upper  boundary: 


w  «=  0  .  uz  *  0 

S  '  «  cons ' t  w/r  to  x  and  t 

n!  =  Ole'  +  ^  [wzz  +  a(H)wz  J 


Lower  boundary: 


Lateral  boundaries: 


w  «  0,  u  «>  0 

G'(x,t)  •>  specified  function 
TT /  *  Cie'  +  [w2Z  +  a(0)wz] 

u  *  0,  w„  “0 

e„'-«  0 

tt'  -  0 


The  conditions  on  tt/  at  the  upper  and  lower  boundaries  follow 
directly  from  8)  in  conjunction  with  the  other  conditions  assumed  at 
these  levels.  In  the  same  way,  certain  other  conditions  follow  as 
results  from  7)  -  10).  These  are  the  following. 


Upper  boundary: 

w2  *■  -ux 

wzz  ■  -a(H)wz  ,  6/z  -  -a(H) 0/ 

Lower  boundary: 

w,  ■  0 

Lateral  boundaries: 

-  0 
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These  conditions  permit  the  formulation  of  all  terms  in  the  various 
equations  of  the  system  in  finite  difference  form.  The  special  treat¬ 
ments  needed  to  obtain  certain  of  the  terms  in  the  balance  equation  will 
be  presented  later. 

The  pressure  computation.  When  the  balance  equation  is  formulated  in 
finite  difference  fora  and  making  use  of  the  derivative  boundary  condi¬ 
tions,  it  is  found  that  the  determinent  of  the  matrix  of  this  elliptic 
equation  is  zero.  As  a  consequence,  the  values  of  Tt'  which  result  from 
solution,  regardless  of  the  procedure  used,  will  be  only  relative  values. 

This  would  constitute  no  particular  problem  if  we  were  concerned 
only  with  dry  convection.  The  pressure  (actually  revalues  are  needed 
only  to  evaluate  the  pressure  force  terms  in  the  equations  of  motion, 
and  since  these  terms  depend  upon  differences  of  pressure,  relative 
values  would  be  acceptable.  However,  as  pointed  out  before,  for  the 
deep,  moist  convection  model  it  will  be  necessary  to  know  unique  values 
of  tt'  at  every  grid  point  in  order  to  determine  the  saturation  mixing 
ratio. 

The  unique  values  are  obtained  simply  by  specifying  the  value  of 
tt'  at  some  point  on  the  boundary  for  all  time.  The  difference  between 
this  value  and  that  computed  for  the  same  point  then  constitutes  a  cor¬ 
rection  factor  to  be  applied  at  all  grid  points. 

Another  problem  related  to  the  pressure  computation  became  apparent 
when  we  were  trying  to  establish  the  initial  pressure  distribution.  The 
basic  form  of  the  balance  equation  may  be  written 

V3™'  +  a(z)nz  «■  C1[0/  +  a(z)e']  +  F(x,z)  12) 
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where  F(x,z)  represents  the  motion  terms  on  the  right-herd  side  of  11). 
We  assume  that  Initially,  the  0'  and  n'surfacea  are  horizontal  and  there 
Is  no  motion.  Thus,  at  the  initial  time,  12)  becomes  simplified  to 

tt/2  +  a(z)TT/  -  Cl[e/  4-  a(s)  0 '  ]  13) 

from  which  one  should  be  able  to  compute  the  initial  tt'  distribution, 
given  the  initial  0'  distribution.  However,  at  the  initial  time  the 
pressure  field  is  hydrostatic  and  thus  one  also  should  be  able  to  com¬ 
pute  the  tt'  field  from 


tt/  -  C10'  14) 

or 

tt/2  -  cie/  15) 

Equation  14)  can  be  solved  by  straightforward  integration,  assuming 
a  value  of  tt'  at  either  the  top  or  the  bottom  of  the  system.  Equation 
15)  can  be  put  in  finite  difference  form  and  solved  either  by  relaxation 
or  as  a  system  of  simultaneous  equations.  Equations  14)  and  15)  give 
the  same  results.  However,  13)  which  simply  represents  a  combination  of 
14)  and  15)  did  not  give  the  same  results  when  solved  as  a  centered- 
dlfference  system.  This  implied  that  once  the  motion  wss  established, 
at  which  time  12)  is  the  only  available  diagnostic  relation  for  the 
computation  of  tt',  we  could  not  expect  to  get  correct  values  of  tt'. 

The  answer  to  this  problem  was  found  by  working  with  various  teat 
cases.  It  was  found  that  in  order  for  13)  to  give  correct  results,  it 
was  necessary  to  express  rrz'7  as  a  centered  difference  at  the  level  in 
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question  while  expressing  tt/  ss  s  difference  over  the  layer  above  or 
below  the  level  of  Interest.  At  the  upper  and  lower  boundaries,  the 
boundary  condition  for  tt/  is  applied  in  the  usual  manner. 

This  procedure,  however,  has  a  disadvantage  in  that  it  introduces 
different  orders  of  truncation  error  in  the  finite  difference  formula¬ 
tion  of  tt/  and  tt/2  .  Thus,  it  was  necessary  to  find  a  way  to  still 
express  tt/  as  a  centered  difference  yet  not  introduce  incorrect  valu 
of  tt'  at  the  Initial  time.  Also,  special  consideration  had  to  be 
given  to  the  formulations  at  the  second  level  in  the  model  due  to  the 
unequal  apace  increments  above  and  below  this  level.  The  manner  in 
which  these  problems  were  handled  will  be  described  in  the  next 
section. 


C.  Numerical  Procedures 

An  extremely  Important  aspect  of  the  modeling  problem  is  the 

choice  of  the  numerical  procedure  which  is  to  be  used  for  the  solution. 

The  procedure  must  be  stable,  it  should  have  as  small  a  truncation  error 

as  possible  and  should  not  be  such  that  it  requires  an  excessive  amount 

of  computer  time  or  exceeds  the  storage  capabilities  of  the  machine. 

Generally,  some  compromise  of  features  will  be  required. 

Numerical  schemes  which  have  been  used  successfully  in  the  past  in 

convection  problems  of  the  sort  considered  here  have  been  the  upstream 

differencing  procedure  used  by  Ogura  (1962,  1963)  and  Orville  (1964, 

1965,  1968),  and  the  two-step  Lax-Wendroff  (1960)  method  used  by 
✓ 

Amason  et.al.  (1968).  Molenkamp  (1968)  has  investigated  the  features 
of  these  schemes  and  others  by  the  application  of  each  to  the  same 
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field  in  a  situation  for  which  the  correct  results  were  already  known. 
His  results  confirm  a  feature  that  we  had  noted  which  is  that  the  up« 
stream  scheme  has  a  built-in  fictitious  diffusion  which  greatly  smooths 
the  field.  On  the  other  hand,  the  lsollnes  of  the  field  were  advected 
at  about  the  proper  rate.  The  two-step  Lax-Wendroff  scheme,  on  the 
other  hand,  conserves  maxima  and  minima  of  the  field  quite  well  but 
advects  the  lsollnes  much  too  Blowly.  Other  schemes  discussed  by 
Molinkamp  have  other  draw-backs  for  grids  of  many  points  because  of 
storage  requirements  or  running  time. 

The  programs  which  we  have  written  up  to  this  time  have  been  based 
on  a  two-step  upstream  difference  scheme  which  will  be  described  below. 
General  features.  We  consider  an  equation  of  the  form 


3F  .ftP  .  „daF 

H  "  mAbi  +  KdZ 


16) 


where  the  first  term  on  the  right-hand  side  represents  advection  and 
the  second  term  represents  diffusion.  In  finite  upstream  difference 
fora  16)  may  be  written  as 


•r  ■  Fi  -  oj  -  rii)  +  <Fi« + Fii  -  2fD  i7> 

where  x  •  nAx,  t  ■  JAt,  r  *  y  -  ■jjl,  A  >  0,  and  Ax,  At  are  the  grid 
and  time  increments,  respectively.  With  simple  manipulation,  17)  can 
be  written 


FJ+1 


- 


K'll  -  'll)  +  Q  +  VIFn 


n+1  +  Fn-1 


18) 
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If,  In  turn,  18)  Is  converted  beck  into  differential-equation  form  by 
using  the  Taylor  Series  through  the  second-order  terms,  one  obtains 


Thus,  it  is  seen  that  representing  the  advection  term  in  16)  as  an  up¬ 
stream  difference  actually  amounts  to  the  use  of  19),  thereby  intro¬ 
ducing  an  additional  diffusion  term  with  the  coefficient, 

'  |A|4*(l--40 

For  values  of  ]a|.  Ax,  and  At  which  are  typical  of  the  convection 
problem,  this  coefficient  is  of  the  same  order  of  magnitude  as  the  real 
eddy  coefficient,  K.  This,  of  course,  presents  a  serious  problem  when 
one  is  interested  in  determining  the  influence  of  real  diffusion  on  con¬ 
vective  development.  It  can  be  shown  that  the  same  influence  arises 
with  the  Lax-Wendroff  Scheme;  however,  from  Molenkamp's  results  th« 
effect  apparently  is  much  smaller.  Presumably,  this  is  because  the 
equation  which  involves  the  ficticious  diffusion  is  used  only  for  every 
other  time  step.  Possibly  this  problem  could  be  overcome  by  use  of  a 
real  eddy  coefficient  whose  magnitude  Increased  as  a  function  of  the 
motion,  i.e.,  along  the  lines  of  the  non-linear  diffusion  discussed  by 
Kuo  (1962).  However,  the  most  satisfactory  answer  would  be  a  numerical 
scheme  which  did  not  introduce  a  spurious  diffusion.  We  have  under 
study  an  implicit  scheme  which  may  prove  to  be  superior  to  any  of  the 
schemes  mentioned  above. 
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The  two-step  upstream  difference  scheme.  In  this  procedure  17)  Is 
computed  In  two  steps  for  each  time  Increment.  First,  intermediate 
values  of  F,  called  F',  are  determined  on  the  basis  only  of  diffusion. 
The  intermediate  values  then  are  used  to  evaluate  thr  advection  term 
and  to  obtain  the  final  predicted  value  of  F.  In  two  dimensions, 
the  system  becomes 


'(FifM 


+  F 


J 

n-l,A 


+  F 


J 

n,jN-l 


+  F; 


n,jt-l 


-  r>[( F'  Xi  •  C F'  Xi.j] 

(  f'  Xi-i] 


20) 


where  x  »  nAx,  z  «  4As,  t  ■  jAt,  Ax  *•  As,  y  •  ri  -  U~,  va  “  Wg|, 
and  U  and  W  are  constant  advection  velocities. 

Let, 


^  ^iknAx  eimXAx 

C  f'  X*  ■  ^  «anix  «lm'“ 


21) 


where  k  and  m  are  wave  numbers  in  the  x  and  z  directions,  respectively 
and  i  -  /-l.  When  the  relations  in  21)  are  substituted  into  20),  one 
finds 


B 

A 


l  -  4-y(^sir?  +  sins 


mAx'N 

2  J 


a  r 
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and  the  stability  criterion  becomes 

1*1  -  r"{ci-(ri+r8)  +  ri  cos  kAx  +  ra  cos  mAx? 

+  [rx  sin  kAx  +  ra  sin  mAx]®J-  <  1 


22) 


It  is  seen  that  the  criterion  Is  a  function  of  wavelength;  however, 
the  result  for  tha  shortest  wavelength  which  can  be  resolved  by  the 
numerical  scheme  Is  found  to  hold  for  all  longer  wavelengths.  Thus,  If 
we  assume  wavelengths  of  2Ax  in  both  the  x  and  z  directions,  the 
condition  becomes 


|*l  -  rni^n+ta)]2  £  l  23) 

where 


r  -  1  -  8Y 

In  the  event  that  the  motion  Is  zero,  which  is  the  Initial  condition  in 
our  caces  unless  an  ambient  wind  field  exists,  ri  and  ra  are  zero  and 
we  must  have  | T j  £  1.  This  condition  is  met  if  y  s  i;.  We  require  that 
this  condition  shall  always  be  met.  Under  these  circumstances,  the 
remainder  of  the  criterion  may  be  stated  as 


|i  -  2(n+ra)|  £  l 

or  24) 

lri  |+|m  |  £  1 

The  absolute  value  signs  have  been  applied  in  24)  to  account  for 
the  modifications  that  arise  in  the  formulations  of  the  advective  terms 
when  the  velocities  are  negative. 
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If  we  let 


M-W-Jlfi'* 

where  V*  la  the  magnitude  of  the  maximum  velocity  vector  In  the  field, 
the  final  results  are  that 

_  Ax8  . 707  Ax 

At  s  and  At  £ - 25) 

A  schematic  graph  of  the  conditions  in  25)  iB  as  follows 

At 

Ax 
4k 


The  allowable  values  of  At  fall  in  the  hatched  area.  It  follows  that 
in  the  initial  stages  of  a  computation.  At  may  be  relatively  large  but 
it  must  be  monitored  and  reduced  as  the  motion  develops  and  increases 
In  magnitude. 

The  truncation  error  for  the  system  given  by  20)  was  determined  in 
the  same  manner  as  that  given  by  Williamson  (1966),  i.e.,  it  is  defined 
as  the  difference  between  the  finite  difference  approximation  equation 
and  the  original  differential  equation.  The  truncation  error  was  found 
to  be 

T  -  0  (At)  +  0(Ax) 
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The  ftntce  difference  relations  for  the  two- dimens tonal  case.  We  con¬ 
sider  here  Che  finite  difference  fonnulstiona  of  7),  8),  10)  end  11). 
Rather  than  setting  up  the  general  relations  for  some  point,  it  will 
be  more  suitable  to  examine  Che  various  terms,  since  in  some  Instances 
rather  special  formulations  are  invol red. 

It  will  be  convenient  to  drop  the  primes  on  tt'  and  in  vhat 

follows  and  to  let  j,k  stand  for  x,s,  where  j  ■  1,2, . je  and  k  - 

1,2, . ke.  The  values  j  ■  1,  Je  and  k  •  1,  ke  denote  points  on  the 

boundary.  The  else  of  the  model  is  easily  varied  by  changing  the 
values  of  Je  and  ke  in  the  computer  program.  It  will  be  recalled  that 
Ax-Aa^lOO  m  except  that  At^lO  m  between  k>l  and  k»2.  In  the  following, 
we  shall  use  Ax  for  a  100-m  increment  and  At  for  the  10-m  increment . 

The  time  derivatives  in  7),  8)  and  10)  are  all  evaluated  as  for¬ 
ward  differences  In  time  in  accordance  with  the  discussion  given  above 
on  numerical  stability.  On  the  other  hand,  the  exact  form  of  the 
advection  terms  depends  upon  the  signs  of  u  and  w.  If  we  let  a  stand 
for  u,  w,  or  0,  the  advection  terms  may  be  written  as  in  the  following 
table. 

Term  u>0,w>0  u<0,w<0 _ 


da 

dx 


,  Ihi LlIizLJs  u  °W :  ZlA 

J  ,k  Ax  J,k  Ax 


da 


w  .1 « k  ~  j  ,  k- 1  k  >  2 
J  ,k  Ax 


j»2  At 


J.Jc 


'  gJ,k 
Ax 
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Because  of  the  boundary  conditions,  all  of  tha  advaction  terms  vanish 
at  tha  comers  of  the  region  and  along  the  bottom  boundary.  At  the  top 
only  uug  is  computed  and  at  the  lateral  boundaries  wwz  and  w  6Z  are 
computed. 

The  Laplacian  terms  in  7) ,  8)  and  10)  are  evaluated  in  the  usual 
manner  at  the  interior  points  as 

-  gl+l.k  +  Vl,k  +  "jjte»l  +  gj,k-l  ~  **J,k  26) 


except  at  the  second  level.  At  this  level  one  obtains  by  use  of  the 
Taylor  Series  expansion  through  the  second  order  terms, 


“l+l.k  +  “l-l.k  "  2“l.k 

-1— • — Er - ^  +DN4o,j  mi 


-  DN^Jfk.i  - 


27) 


WheM  "  toiCM&V  DN5  “  da  (Ax+Az)  *  and  DN6  ‘ 

Expressions  then  are  needed  for  uz  z  at  the  upper  boundary  and  wx  x  and 
8XX  at  the  lateral  boundaries.  The  boundary  conditions  described  in 
Section  B  permit  these  derivatives  to  be  expressed  as 


2(ati  -  tto) 
Ax2 


28) 


where  n  is  in  the  normal  direction  to  the  boundary,  oi  is  the  value  of 
a  at  the  point  just  Interior  to  the  boundary,  and  a0  1®  the  boundary 


value . 


The  flrsC  derivative  terse  In  the  diffusion  quantities  of  7),  *) 


end  10)  ere  treated  aa  centered  differences  and  take  the  font 


except  at  level  two  where  again  a  special  formulation  is  required. 
Using  the  Taylor  Series  expansion  through  the  second  order  term,  one 
obtains  for  this  level  the  expression 


.el  -  «• 


J,k+1  *  “^j.k-l  +  “^J.k 


where  DNl  ■  -r — 7^~"?"*r,  DH2  «  .  ■  — r ,  and  .‘dJ3  • 

Ax(Ajn-As)  &s<A*tA*) 

The  pressure  force  tens  in  7)  is  spproxlestsd  by 


at  all  interior  points  and  along  the  upper  boundary.  The  pressure  tens 
In  8)  also  is  represented  ay  s  centered  difference;  however,  this  choice 
must  be  considered  in  conjunction  with  the  wey  the  buoyancy  tern  In 
formulated  and  also  must  oe  consistent  with  the  way  the  teem  are 
handled  in  11).  Thus,  we  shall  take  this  up  in  the  following  discussion 
of  11) . 

The  pressure  computation  must  be  carried  out  at  the  boundary  points 
as  well  as  at  the  interior  points  of  the  grid.  Therefore,  a  mans  must 
be  found  to  express  all  of  the  teres  In  11)  which  involve  a  first  or 
second  derivative  that  Is  taken  normal  to  a  boundary  and  which  are  not 
already  accounted  for  in  the  boundary  conditions  stated  In  Section  5. 
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Specifically,  at  the  lateral  boundaries  ve  need  ttxx  and  u* .  Using 
the  conditions  at  these  boundaries  that  tt  ■  0  and  u  ■  0,  we  can  express 
r0  in  the  same  form  as  that  shown  in  28)  while  ux  becomes 


a.  ±  a.,:.. 

5x  Ax 


U£ 


32) 


where  the  +  sign  applies  at  the  left-hand  boundary  and  the  -  sign  ap¬ 
plies  at  the  right-hand  boundary. 


At  both  the  upper  and  lower  boundaries  we  must  have 


2  2 


These 


expressions  are  readily  obtained  in  centered  difference  form  by  using 

the  upper  boundary  condition  tt2  ■  C10,  and  the  lower  boundary  condition 
Y  2w 

ttz  »  C16  +  -jpg)  Accordingly,  we  obtain  at  the  upper  boundary. 


and  at  the  lower  boundary, 


2AxC19. 


ke 


33) 


^  -  <v ;  "i,0 ,  ^ 


w 


T~  Cclei  l 
Az  V  J.l 


34) 


Statements  for  02  also  are  required  at  the  upper  and  lower  boundaries. 
An  attempt  to  use  the  thermal  equation  10)  as  a  diagnostic  relation  to 
obtain  6Z  at  these  boundaries  failed  because  the  resulting  balance 
equation  could  not  provide  a  hydrostatic  distribution  of  pressure  at 
the  initial  time.  Consequently,  these  derivatives  are  approximated  by 
inward  differences.  Then  at  the  upper  boundary  we  have 

OS 
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and  at  the  lower  boundary 


bz 


36) 


The  expression  for  6Z  at  level  k  ■  2  had  to  be  derived  by  special  means. 
We  consider  the  region  at  the  Initial  time  when  there  is  no  motion,  the 
0  surfaces  are  horieontal,  and  the  heating  at  the  bottom  is  sec  equal  to 
zero.  At  this  time,  the  pressure  Is  hydrostatic  and  15)  could  be  used 
to  compute  the  distribution.  An  integration  of  15)  with  respect  to  z 


yields 


■  <eke  -  0 

which,  of  course,  agrees  with  the  boundary  conditions.  If  the  integral 
of  02  is  expanded  using  the  trapezoidal  rule,  we  must  have 


0ke  '  91  “  L^1  +  (0^)8]t 


+  f%«z>3  +  <ez>3  +••••+%  <Mke]Ax 


On  the  right-hand  side  we  use  35)  and  36),  and 


k  ^  1,2  or  ke 
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Thus,  37)  is  an  equation  for  (8 t)a •  After  cancellation  of  many  terms 
and  rearrangement,  the  result  is 


('M'\  ^  ea  -  &L 
Ax  +  Lz 


38) 


The  final  quantities  to  be  considered  in  11)  [and  also  in  8)]  are 
the  terms  a(z)TTr  and  Cla(z)0.  These  terms  are  related  to  one  another 
at  the  upper  and  lower  boundaries,  as  already  noted.  Elsewhere,  as 
pointed  out  in  Section  B,  it  was  found  that  at  the  initial  time  the  only 
way  the  balance  equation  11)  could  give  the  same  results  as  the  hydro¬ 
static  equation  14)  was  if  rr2  and  8'  were  evaluated  for  the  layer  above 
or  below  the  level  of  interest.  Let  the  level  of  interest  be  denoted  by 
k.  Then  the  hydrostatic  equation  gives 


and 


\+l  ~  nk  .  (9k+l  +  9kl 
Ax  “  Cl  2 


^  ~  Vl  .  C1  6k  +  ?fel| 

Ax  C1  2 


39) 


Thus,  one  can  satisfy  the  hydrostatic  relation  at  level  k  by  combining 
these  two  equations  into  the  single  relationship 


At  level  two,  thla  becomes 


In  this  manner,  nz  is  expressed  as  a  centered  difference  while  9  Is 
represented  by  a  weighted  average  in  8)  and  11).  After  the  initial 
moment,  the  heating  at  the  ground  disturbs  the  hydrostatic  distribution 
of  pressure  which  in  turn  gives  rise  to  accelerations  in  7)  and  8). 

Equation  11)  is  solved  for  tt  by  an  iterative  relaxation  technique 
(Llebman  process).  The  finite-difference  form  of  11)  at  Interior  points 
(3«k£ke-l)  is 

A?  Gj+l.k  +  nj-l,k  +  nj,k+l  +  "j.k-l  '  4TTj,k) 

42) 

+  ~2tot  0rj  ,k+l  ‘  ^j.k-O  ■  Fj,k 

All  of  the  terms  on  the  right-hand  side  of  11)  are  included  in  F  , 

J 

(x,r).  At  the  boundaries,  the  special  formulations  of  derivatives  given 
above  must  be  used. 

At  the  level  k“2,  11)  takes  the  form: 


Etcher  42)  or  43)  can  be  solved  for  n  .  ,  so  thee 

"j.k  •  "j.k 

where  H  contains  all  ocher  terns  in  42)  or  43) . 

The  next  approximation  for  rr  la  found  as  tt*: 

<"j.k  •") 

Here,  m  is  the  underrelaxatlon  factor  which  gives  convergence  in  the  case 
0  <  m  <  1.  The  convergence  of  the  described  relaxation  process  was 
established  by  several  numerical  experiments  with  various  values  of  m. 

The  fastest  convergence  was  found  to  occur  for  m  •  0.5. 

Integral  Constraints.  It  is  standard  procedure  in  modeling  work  to 
form  the  mechanical  energy  equation  for  the  system.  For  a  closed  system, 

such  as  is  usually  dealt  with,  this  amounts  to  a  statement  of  the  con¬ 

servation  of  the  sum  of  the  kinetic  end  available  potential  energies. 

The  relation  may  then  be  used  as  a  check  on  the  accuracy  of  the  numerical 
integration. 

In  the  present  study,  the  system  is  open  to  energy  input  which  is 
reflected  in  a  somewhat  more  complicated  relation  for  the  energy  equation. 
The  derivation  follows  in  the  usual  manner  from  the  equations  of  motion. 
For  dry,  deep  convection,  v’e  obtain  from  1)  -  3)  the  relation, 

i(^).  -q.8C.VTT'  +  wgf  +  K*  [u^u+v^vH^v] 

44) 

.  .  f  du  .  5v  ,  dw~i 

+  Km  a(z)  [ujj  +  v^  +  v^J 
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When  44)  la  multiplied  through  by  Po  and  the  continuity  equation  is 
used,  ve  obtain 


■^PoK  -  -  -^op0yc  -Cp  ®7#Pott '<£  +  |  7.po0/zC 

-  +  PoK*  a ( z)  [ u^3 ufvV3  v+w'7® w] 


45) 


,  V  ,  v  5K 
+  PoK*  a(z)  ^ 

ua  -t-v^ 

where  K  - r~~ —  •  The  fourth  term  on  the  right  may  be  reexpressed  by 

using  5),  When  45)  is  integrated  over  the  volume  of  the  model  under 
conditions  of  a  zero  ambient  wind  field  and  use  is  made  of  the  boundary 
conditions,  the  resulting  equation  is, 


where  H  stands  for  the  upper  boundary  of  the  system  and  dt>  »  dxdydz. 

The  left-hund  side  of  46)  expresses  the  change  of  kinetic  and  available 
potential  energy.  The  dissipation  terms  on  the  right-hand  side  appear 
because  the  thermal  energy  equation  5)  does  not  contain  a  term  accounting 
for  the  frictional  transformation  of  kinetic  into  heat  energy. 

A  somewhat  simpler  relationship  may  be  derived  from  the  thermal 
energy  equation.  In  this  case,  we  multiply  5)  by  p0  and  integrate  over 
the  volume  of  the  model.  From  this  ve  obtain  after  application  of  the 
boundary  conditions 
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£t  J  pt'e'du  ■  J  (f50^  dxdy  -  J  (pc**  ***y  *7) 


Thus,  In  this  case,  the  change  in  sensible  heat  of  the  region  is  seen 
to  depend  upon  the  hast  fluxes  at  top  and  bottom  of  the  system. 

The  two-dimensional  versions  of  46)  and  47)  follow  directly  by 
setting  v  and  ^  equal  to  zero. 


D.  The  Ambient  Flow  Model 

General  features  This  model  follows  in  a  straightforward  manner 
from  the  anelastic  system  given  previously  by  1)  -  6)  if  the  coriolis 
force  is  incorporated  for  purposes  of  modeling  the  large-scale  flow. 

It  is  assumed  here  that  the  atmosphere  consists  of  a  stationary, 
hydrostatic  state  as  represented  by  the  zero-order  equations  given  by 
Ogura  and  Phillips  (op.cit.).  Superimposed  upon  this  is  a  large-scale 
basic  or  ambient  flow  which  for  simplicity  may  be  taken  as  zonal  and 
steady-state.  This  flow  may  be  modeled  to  include  a  friction  layer  in 
which  vertical  eddy  momentum  transfer  occurs.  Geostrophic  flow  may  be 
assumed  above  the  friction  layer,  and  both  horizontal  and  vertical  shear 
may  be  permitted.  A  convective  perturbation  may  then  be  added  which  is 
ageostrophic  and  of  the  same  order  of  magnitude  as  the  ambient  flow. 

An  x-z  plane  representation  of  the  model  is  given  in  Fig.  II-2. 

It  is  assumed  that  initially  tt*  ■  0,  9/  ■  0  and  only  the  ambient  flow 


exists. 


I 
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Fig.  II-2  Ambient  Flow  Model  (x-y  plane) 

As  in  the  previous  model,  differential  heating  takes  place  along  the 
ground  leading  to  the  development  of  hot  spots.  As  shown  in  Fig.  11*2, 
it  is  assumed  that  conditions  in  the  model  section  (solid  boundaries)  are 
repeated  in  the  sections  upstream  and  downstream  relative  to  the  basic 
flow. 

System  of  equations.  In  the  following  discussion  we  let  the  subscript 
1  refer  to  ambient  flow  and  the  subscript  2  refer  to  the  convective 
fields.  The  dependent  variables  of  the  system  may  then  be  represented 


u  -  ux  (y,z)  +  mb  <x,y,z,t) 

V  -  Vs  (x.y.z.t) 

w  -  vs  (x,y,z, t) 

it'-  tti  (y,s)  +  tt2  (x,y,z,t) 

e'-  0i  (y)  +  02  (x,y,z»t) 


It  may  be  noted  that  the  lapse  rate  for  the  ambient  flow  is  adiabatic 
but  there  is  a  north-south  temperature  gradient.  For  a  two-dimensional 


version  of  this  model,  0i ,  may  be  set  equal  to  2ero  which  means  that  the 
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ambient  flow  i8entropa  16  given  by  9,  and  ttj  ■  Tfe.  From  1)  *  5),  the 
ambient  flow  equations  are  given  by: 


0  "  0  ^  -fui  +  Kh  +  a<*>  fe] 

0-.Cp®^+8^  49> 


0 


0 


When  the  statements  in  48)  are  substituted  into  1)  -  6)  and  cancel 
lations  are  made  in  accordance  with  49),  we  obtain  for  the  three- 
dimensional  model: 

+  +  «  ft  ■  ♦  ■">  fc] 

50) 

*fc  *  T&  +  n  ^  '  *■[’’’ «>  +  •<*>  W] 


?°Po"S  +  *  0 


+  a(nB)2  ■  Cl[(Qa)j  +  *0a  ]  +  “g  [  "IT  tt^  +  IT  J 

+  CllC (v% ) x  (us-ui)2  -  (ua),  (ws):  +  (va)x  (ua-ui)Y 

-  (ua)x  (va)y  +  (va)2  (Wa)y  -  (va)Y  (w&)2  ]  50) 

These  equations  with  appropriate  assumptions  may  fojm  the  basis  for 
several  different  dry  models. 

1)  A  tonal  basic  flow  with  vertical  shear  (two* dimensional) 

2)  A  ronal  basic  flow  with  vertical  shear  (three-dimensional) 

3)  A  zonal  basic  flow  with  horizontal  and  vertical  shear 

(three-dimensional) 

It  is  possible,  also,  with  some  alterations  to  Include  the  influence  of 
horizontal  convergence  in  the  ambient  flow.  The  possibility  of  moist 
convection  will  be  considered  at  a  later  point. 

Boundary  conditions.  The  boundary  conditions  for  this  model  are  in  moat 
respects  the  same  as  those  for  the  model  with  zero  ambient  flow.  The 
primary  difference  occurs  In  the  conditions  in  the  x  direction.  As 
indicated  In  Fig.  II-2,  all  dependent  variables  are  repeated  periodically 
in  the  x  direction.  Thus,  Q(x)  -  Q(jm-L),  where  Q  is  any  one  of  the 
variables.  In  addition,  the  following  conditions  are  assumed. 

Upper  boundary:  Vs  -  0,  (u3>j  •  0,  (uj.)z  "0,  (va)j-  0 

0e  ■  cons' t.  w/r  to  x,  y,  and  t, 

(tt2),  -  Clfe  +  fg  (ws)22  +  (H)  (W3  )j  j 
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Lower  boundary: 


Wb  ■  0,  ua  »  0,  ui  ■  0,  va  ■  0 
9a  (x,y,t)  ■  specified  function 
(na)j  ■  Cl9a  +-^%r(wa)2J  +  e(0)(wa),1 

North-South  boundaries:  (ws)y  ■  0,  (ua)Y  ■  0,  va  ■  0, 

(ux)j  ■  specified  function 

(9a)v  »  0 

(rra)y  ••  0 

When  these  boundary  conditions  are  applied  to  50),  the  following 

additional  relations  are  found; 

Upper  boundary: 

Lower  boundary: 

North  South  boundaries: 

E.  Molat  Convection 

General  Considerations.  This  model,  especially  with  the  inclusion 
of  an  ambient  wind  field,  is  the  ultimate  goal  in  our  work  although  it 
has  received  the  least  attention  in  this  early  phase  of  the  project. 

Dr.  P.  Das  has  concentrated  most  of  his  effort  on  the  development 
of  relations  wh.'ch  will  incorporate  Che  microphysics  of  the  phase  changes 
in  a  convective  cloud.  A  summary  of  this  work  will  be  given  at  a  later 
point  in  this  section.  However,  for  the  time  being  we  again  turn  to  the 
anelastic  system  devised  by  Ogura  and  Phillips  (op.cit.),  which  is  suit¬ 
able  for  modeling  the  macroscale  aspects  of  cumulus  development,  assuming 
a  reversible,  saturation  adiabatic  process. 


(U»)x  +  (V3>y  +  (ws);  -  0 
w,  t  *  -a(H)w; 

(Wa)2  -  0 

(vs  )  y  y  “  0 


Tha  applications  of  this  aystsai  (Ogura,  1963;  Orville,  1965  and  1966) 
hava  always  boon  to  oballow  convactlon,  apparently  to  ovoid  tho  coop  11  ca¬ 
tion  which  tirlaaa  In  tho  computation  of  tho  saturation  nixing  ratio  duo 
to  ita  dependsnca  upon  proaauro.  As  aXow.i  hy  Ogura  (1963),  an  approx¬ 
imation  to  tho  aaturacion  vapor  prasaura  la  g,lve.i  by 

-  ••<*>[»  +  i^r  ♦  *G^y7  T"] 

whore  L  la  tho  latont  hoot,  By  la  tho  apocific  gaa  conatant  for  vapor, 
and  T0  la  tho  roforonce  atmosphere  temperature.  For  shallow  convection, 
t'  Is  the  deviation  of  temperature  from  the  roforonce  value;  however, 
for  deep  convection  T'  la  related  to  the  potential  temperature  and 
pressure  deviations  through  the  relationship, 

t'  e'  n' 

T0  “  0  no 

More  will  be  said  about  this  shortly. 

System  of  equations.  In  this  system  the  predictive  equation  for 
potential  temperature  deviation  Is  replaced  by  one  for  the  specific 
entropy  of  the  moist  fir.  Also,  a  term  is  added  to  the  third  equation  of 
motion  which  accounts  for  the  negative  buoyancy  of  the  liquid  portion  of 
the  water  substance.  The  complete  system  is: 

•  -CpOn*'  -r  Kn^u  +  a(z)\iz]  51) 


36 


52) 


g  -  -C,&nr*  +  H,[^v  +  a(z)v2] 

g  -  -q,  ®rr2'  +  ^  -  g  <L  +  KmC^w  +  a(z)w2  ]  53) 

V°p0$  ■  0  54) 


g  -  K„[V2q  +  a(z)qz  ]  where  q  •  q*  +  qL 


55) 


g  ■  Kp[Va<p  +  a(z)epj  ] 


flil 


e.' 


cp  -  -*#*■  +  <r-^  q*  or  cp  =  0S  + 


Q.T, 
0'  IT-N 


q»s 


[>♦  *X  (®'+  O  +  <  vc) 
(!’♦  0  ] 


56) 

57) 


58) 


where  qv  and  qL  are  the  mixing  ratios  for  vapor  and  liquid,  cp  Is  the 
specific  entropy  of  the  moist  air  (the  entropy  of  the  liquid  is  neglibly 
small),  and  9jNJ  and  0,  are  the  potential  temperature  deviations  for  un¬ 
saturated  and  saturated  air,  reaper  ively.  Equation  58)  expresses  the 
saturation  mixing  ratio  where  p  «  Prr1^  »  P(tt0  +  n')1^  .  The  divergence 
or  balance  equation  for  the  system  is  given  by, 


V=n'  +  a(z)n/  -  Cl[9/  +  a(s)0']  -  C3[(qL)2  +  a(z)q,J 


TT0 


r-c7  w 
L  4  T?o 


59) 


+  Cll[w,u2  -u,wz+vxu,  -u*vY+VjWv  -w2  vY  ] 
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In  order  to  determine  whether  or  not  condensations  has  occured  at 
any  given  grid  point,  one  solves  55)  and  56)  for  a  q  and  <p  at  the  point. 
These  values  are  then  used  in  the  first  expression  of  57)  to  compute 
0jNa  *  The  value  of  tp  again  is  used  in  the  second  expression  of  57) 
which  is  solved  simultaneously  with  58)  to  get  8,'  and  qv!.  If 
8j  <  OJnj,  then  qL  •  0  and  0'  -  0jNS .  However,  if  8*  >  ©Jns ,  then 
0*  ■  0j  and  q^  ■  q  -  qv  s  ■  The  special  problem  that  arises  in  this  pro¬ 
cedure  is  because  of  the  dependency  of  58)  on  rr',  which  must  be  known. 
However,  v‘  is  computed  from  59)  which  depends  upon  d‘  and  q,^  .  Thus, 
a  dilemma  is  created  as  to  how  to  start  the  computation. 

Tentatively,  at  least,  we  propose  to  use  an  iterative  scheme  for 
the  solution,  as  outlined  below.  The  first  step  is  to  assume  for  the 
moment  that  the  air  is  saturated.  Tien  from  the  second  expression  in  57) 


0'  -  0,' 


L@ 
Q>  T0 


s 


Also, 


60) 


<b  s 


If  the  expressions  in  60)  are  substituted  into  59),  the  latter 
becomes, 

^tt'  +  a(z)n'  -  C3-|(jpz  +  acp)  -  Qu  +  aq^)  -  (~~  "  0 

[(qv  sy  +  a<h s  ]  +  qvs}  61> 

,  i  r-c7  w 

L  4  Ti0 


+  T»I  ■  C9W2 
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+  Cl  1  [  (wxuz  -uxw2  )  +  (vxUy-UjjVy)  +  (vz  wY  -w2  vy )  ]  61) 

Equation  61)  then  expresses  the  pressure  field  for  saturated  air;  how¬ 
ever,  it  holds  equally  well  for  unsaturated  air  if  qYS  is  replaced  by 

q- 

Now,  at  the  end  of  a  given  time  step  we  have  the  fields  of  q  and  cp 
but  we  do  not  yet  have  9',  v'  and  q,_ .  Assume  the  old  value  of  tt'  as  a 
first  estimate  in  58)  and  solve  58)  simultaneously  with  the  second 
expression  in  57)  to  obtain  a  first  estimate  of  the  98  and  qv s  fields. 
Substitute  the  values  of  cp,  q  and  qvs  into  61)  to  get  a  new  estimate  of 
the  tt'  field.  Use  the  new  estimates  of  tt'  in  58)  and  solve  as  before 
for  second  estimates  of  9,'  and  qv .  This  procedure  may  be  repeated 
until  the  desired  degree  of  accuracy  is  obtained  for  tt',  0s'  and  qY8, 
assuming  that  convergence  occurs. 

Now  compute  9^3  from  the  first  expression  in  57)  using  cp  and  q, 
and  compare  8^3  with  9S'.  If  0,'  >  0jN5 ,  then  saturation  has  occurred 
and  the  last  values  of  tt  ,  6/  and  qY8  are  the  desired  values,  o^  can 
be  obtained  explicitly  from  qt  ■  q  -  qY8.  If  ®s  <  ®Jns»  t'len  saturation 
has  not  r  curred  and  08  is  the  desired  value  of  9'.  tt'  is  then  computed 
from  61)  by  setting  qv ,  »  q. 

It  is  evident  that  this  model  will  require  considerable  storage 
capacity  in  the  computer  and  the  procedure  will  involve  large  amounts 
of  computer  time  per  forecast  period.  This  would  be  true,  particularly, 
if  the  model  were  modified  to  include  a  shearing  ambient  wind  field  with 
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precipitation  and  in  three  dimensions.  Therefore,  it  well  nay  be  that 
we  shall  have  to  restrict  our  experiments  to  relatively  simple  cases 
when  we  get  to  this  model. 

■T.ll?.P^r°P.hy.Bica  of  moist  convection.  At  this  point,  we  give  a  brief 
reviev  of  Das's  work  which  is  concerned  with  certain  aspects  of  phase 
change  and  liquid  water  storage  in  a  convective  cloud. 

A  theoretical  and  numerical  study  of  moist  convection  is  complicated 
by  the  various  physical  factors  related  to  the  nature  of  phase  transition 
of  water  substances.  In  the  classical  approach  to  the  moist-adiabatic 
process  one  assumes  that  all  the  water  vapor  above  the  saturation  value 
condenses  and  that  in  the  presence  of  liquid-water  particulates  the  cloud 
air  is  always  just  saturated.  Even  this  simple  assumption  entails  serious 
difficulties  in  numerical  modeling  although  unfortunately,  it  can  afford 
only  a  very  simple  picture  of  moist  convection. 

In  order  to  examine -what  is  involved  in  a  more  tv...  ,.c  modeling  of 

moist  convection,  a  reexamination  of  the  thermodynamic  equation  was  taken 
up  for  an  alternative  form  which  would  "ake  care  of  the  microphysical 
implications  of  the  process  of  condensation  and  evaporation  in  convective 
drafts.  The  alternative  picture  arrived  at  replaces  the  moist-adiabatic 
equation  by  a  set  of  three  equations,  viz.,  1)  a  moisture  equation. 


dr 


2)  th*  thermodynamic  equation, 
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d<t>d  ,  L  dgu 
dt  T  dt 

and  3)  an  equation  of  continuity  for  N^, 


dNr 

dt 


—  V.  (n  )  -4  Cn  %?) 
p,  r  r  o  v  r  dt„ 
rd  r 


which  Includes  the  effect  of  condensation-evaporatior  but  not  coagulation. 
The  symbols  have  the  following  meanings: 


d  B  “*  -* 

-tt  t—  +  V>V  where  V  is  the  velocity  of  air  motions, 
<1  c  ■  o  t 


V  ■■  velocity  of  the  particulate  products  of  condensation  relative  to 
air  motions, 

§v  “  humidity  mixing  ratio, 

N  ■  number  of  particles  taking  part  in  the  processes  of  condensation 
and  evaporation  per  unit  mass  of  dry  air, 

in  ■  effective  water  mass  of  the  particles  of  condensation  of  radius  r, 

cp^  ■  entropy  per  unit  mass  of  dry  air, 

L  ■  latent  heat  of  condensation, 

T  ■  temperature  of  cloud  air, 

Pj  *  density  of  dry  air 


n  ■  p  .N 
r  rd  r 


number  of  particles  of  condensation  per  unit  volume  of  cloud. 


The  equation  for  N^,  which  does  not  include  the  effect  of  coagulation, 
can  be  made  to  do  so  by  including  an  additional  set  of  terms.  Fortunately, 
over  the  last  few  years  this  aspect  of  cloud  processes  has  received  able 
consideration  from  Golovin  (1963),  Twomey  (1964,1966),  Berry  (1965), 
Bartlett  (1966)  and  Warshaw  (1967)  and  the  cloud  dynamicist  will  find  the 
treatment  ready  for  inclusion  in  a  model  when  he  has  one  tractable  enough 
within  the  existing  computation  abilities. 

The  immediate  achievement  of  the  new  thermodynamic  equation  is  in¬ 
corporating  in  one  equation  the  effects  discussed  by  McDonald  (1962), 
Kamburova  and  Ludlam  (1966),  and  Das  and  Subba  Rao  (1968)*.  McDonald 
considers  the  anomalous  thermodynamic  consequences  of  supersaturation 
in  an  updraft  containing  a  population  of  condensation  nuclei.  There  are 
circumstances  available  when  the  adiabatic  ascent  of  a  parcel  can  also 
be  isothermal.  On  the  other  hand  Kamburova  and  Ludlam,  and  Das  and  Subba 
Rao  have  shown  that  a  thunderstorm  downdraft  cc  itaining  all  its  water  in 
the  form  of  drops  will  be  unsaturated  and  its  temperature  lapse  rate 
closer  to  the  dry  adiabatic  when  the  water-content  is  divided  into 
larger  drops.  In  other  words,  the  thermodynamic  system  envisaged  by  the 
new  thermodynamic  adi*  xc  process  dispenses  with  the  saturation 
assumption  of  the  classical  moist  adiabatic  process. 

A  paper  incorporating  the  results  of  the  above  studies  was  communi¬ 
cated  to  the  Journal  of  the  Atmospheric  Sciences.  It  has  been  revised 
in  view  of  the  referees'  comments  and  resubmitted  for  publication. 

♦Unpublished  work 
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Introducing  the  tnlcrophysics  of  cloud  formation  into  the  problem 
will  involve  obvious  complications  but  in  its  simplest  aspects  might 
provide  an  important  consequence.  Since  sub-  or  super-saturation  will 
be  a  matter  of  rule,  the  calculations  will  be  less  sensitive  to  the 
exact  point  where  saturation  occurs. 

A  simple  way  to  start  introducing  the  effects  of  the  microphysics 
is  by  assuming  nuclei  of  a  single  size  with  an  assumed  but  realistic 
Concentration.  This  will  not  be  much  more  difficult  than  working  with 
the  pseudo-adiabatic  form  of  the  thermodynamic  equation.  However,  the 
introduction  of  a  population  of  nuclei,  and  then,  of  coagulation  effects 
will  involve  considerably  more  sophistication  and  computer  time. 

The  details  of  the  application  of  these  concepts  to  the  moist 
convection  problem  are  not  as  yet  worked  out. 

F.  Summary 

The  foregoing  summarizes  the  various  models  which  are  to  be  studied 
and  the  lines  along  which  our  work  will  proceed  in  the  next  two  years. 
The  programs  for  the  dry  convection  model  which  have  been  run  on  the 
computer  to  date  have  been  for  testing  purposes;  therefore,  there  are 
no  computational  results  to  present  at  thi3  time. 

Two  graduate  students  presently  are  doing  their  thesis  research  in 


connection  with  this  project.  Captain  R.  L.  Croft  (USAF)  is  working 
with  the  ambient  wind  field  model  with  vertical  shear.  Captain  J.  L. 
Conley  (USAF)  has  initiated  a  theoretical  study  along  the  lines  of  the 


work  done  by  Kuo  (op.clt.),  Saltzman  (1962),  and  others.  The  latter. 

If  successfully  concluded  for  the  boundary  conditions  with  which  we  are 
working,  should  compliment  the  numerical  computations  and  greatly  enhance 
our  understanding  of  the  results. 
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GLOSSARY 


A 

a(z) 

C 

es 

IF 

g 

H 

A  A  A 

i,  j  ,  k 

1 

J 

k 


a  constant  east-west  component  of  wind 

percentage  rate  of  change  of  reference- 
atmosphere  density  with  height 

three-dimensional  velocity  vector 

specific  heat  of  air  at  constant  pressure 

saturation  vapor  pressure 

three-dimensional  friction  vector 

acceleration  of  gravity 

height  of  top  boundary  of  system 

unit  vectors  in  the  x,y,z  directions, 
respectively 


/-I 

grid  point  number  along  time  coordinate 
or  x  coordinate,  as  indicated  in  text 

grid  point  number  along  z  coordinate 
or  a  wave  number  in  x  direction,  as 
indicated  in  text 


(cm/sec) 

(1/cm) 

(cm/sec) 
(ergs/gm  deg) 
(dynes/ cm2  ) 
(dynes/gm) 

(cm/ sec2) 

(cm) 

(dimensionless) 

(dimensionless) 

(dimensionless) 

(dimensionless) 

or 

(1/cm) 


Kh.K,  .Kp.KtB 


K 

Jt 

L 

m 


exchange  coefficient  for  momentum,  heat, 
water  substance,  and  entropy,  respectively 

a  general  exchange  coefficient 

grid  point  number  along  z  coordinate 

latent  heat  of  vapor'  'tion  for  water 

wave  number  in  the  z  direction 


(cm2 /sec) 

(cm2 /sec) 

(dimensionless) 

(ergs/gm) 

(1/cm) 


m 


r 


effective  water  mass  of  particles  of 
condensation  of  radius,  r  (g™) 
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grid  point:  number  along  x  coordinate 

number  of  condensation  particles  per 
unit  volume  of  cloud 

number  of  particles  Involved  in  evaporation 
and  condensation  per  unit  mass  of  dry  air 

pressure  of  air 

a  reference  pressure  for  air 

mixing  ratio  for  water  vapor 

mixing  ratio  for  liquid  water 

mixing  ratio  for  water  substance 

mixing  ratio  for  water  vapor  at  saturation 

specific  gas  constant  for  dry  air 

specific  gas  constant  for  water  vapor 

the  ratio  UAt/Ax 

the  ratio  WAt/Ax 

time 

absolute  temperature 

reference- atmosphere  absolute  temperature 

deviation  of  temperature  from  the  reference 
value 

component  of  velocity  in  the  x  direction 

a  constant  component  of  velocity  in  the 
x  direction 

component  of  velocity  in  the  y  direction 
che  maximum  velocity  in  the  x,z  plane 


(dimensionless) 

(1/cm3) 

(1/gm) 

(dynes/cm3) 

(uynes/cn? ) 

(dimensionless) 

(dimensionless) 

(dimensionless) 

(dimensionless) 

(ergs/gm  deg) 

(ergs/gm  deg) 

(dimensionless) 

(dimensionless) 

(sec) 

(deg) 

(deg) 

(deg) 

(cm/sec) 

(cm/ sec) 
(cm/sec) 
(cm/sec) 
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component  of  velocity  in  the  z  direction 

(cm/ sec) 

a  constant  component  of  velocity  in  the 
e  direction 

(cm/ sec) 

west-east  coordinate,  positive  eastward 

(cm) 

north-south  coordinate,  positive  northward 

(cm) 

vertical  coordinate,  positive  upward 

(cm) 

a  dummy  variable  for  u,  w,  or  S 

(cm/sec) 
or  (deg) 

grid  increment  of  100  m 

(m) 

grid  increment  of  10  tn 

(m) 

time  increment 

(sec) 

potential  temperature 

(deg) 

reference  -atmosphere  potential  temperature 

(deg) 

potential  temperature  deviation 

(deg) 

potential  temperature  deviation  at 
saturation 

(deg) 

potential  temperature  deviation  for 
unsaturated  conditions 

(deg) 

ratio  of  R/Cp 

(dimensionless) 

volume 

(cm3) 

mixing  ratio  for  water  vapor 

(dimensionless) 

a  pressure  variable,  (p/P)K 

(dimensionless) 

reference-atmosphere  pressure  variable 

(dimensionless) 

pressure  deviation 

(dimensionless) 

reference-atmosphere  density 

(gm/ cm3) 
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density  of  dry  air 
specific  entropy  of  moist  air 
specific  entropy  of  dry  air 


(gm/cm3) 
(dimenalonlaM) 
(cals/gm  dag) 
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III.  ANALYSIS  AND  PREDICTION  OF  ATMOSPHERIC  STRUCTURE  THROUGH  USE  OP 
DUAL  FREQUENCY  RADAR 

V.  E.  Moyer 


A.  Introduction 

The  objectives  of  this  phase  of  the  Themis  program  (Sub-task  b) 


are: 


1)  Observe,  record,  and  analyze  precipitation  echo  strength 
and  patterns  using  two  carefully  calibrated  radars 
which,  while  operating  at  different  wave  lengths,  are 
matched  with  respect  to  pulse-repetition  rate,  pulBe 
length,  receiver  band-pass  characteristics,  antenna- 
beam  pattern,  and  antenna  height; 

2)  Observe  precipitation  intensity  and  amounts  using 
calibrated  raingages  over  a  watershed; 

3)  Observe  accompanying  meteorological  and  electromagnetic 
propagation  conditions  at  and  near  the  radar  antennas 
and  the  watershed  area; 

A)  Investigate  cell  mocion  versus  line  motion; 

5)  Develop  prediction  techniques  for  cell  and  line  motion; 
and 

6)  Construct  conceptual  and  quantitative  models  of 
precipitation  phenomena  and  develop  and  test  rea  istic 
and  effective  prediction  techniques  for  preclpit'  .ion 
parameters  including  rainfall  intensity,  liquid-water 
content,  drop-size  distribution,  surface  visibility, 
and  precipitation  amount, 

B.  Procedures 


Through  the  first  year,  the  work  has  proceded  with  the  above 
goals  in  view.  From  the  initiation  of  the  project,  a  24-hour  radar 
watch  was  set  up  to  record  all  precipitation  echoes  that  appeared 
inside  the  125-mile  and  outside  the  25-mile  range  markers.  The 
observing  technique  that  was  used  was  that  of  sequential 
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gain-reduction  of  the  receivers.  At  the  same  time,  efforts  to  optimise 
the  equipment  went  forward. 

By  the  end  of  August  1968,  the  following  characteristics  of  the 
equipment  had  been  achiw.ed. 


Table  1II-1.  Characteristics  of  Texas 

A&M  Dual-frequency 

Radar  System 

Peak  power,  PT 

AN/CPS-9 

2X105  w 

WSR/TAM-1 

2X105  w 

Pulse  repetition  frequency,  PRF 

186  sec-1 

186  sec-1 

Wavelength,  X 

3.2  cm 

10.3  cm 

Pulse  length,  t 

5.0  psec 

5.0  psec 

Pulse  length,  h 

1.5X105  cm 

1.5X10 5  cm 

Antenna  shape 

Paraboloidal 

Paraboloidal 

Antenna  diameter 

1.52X102  cm 

4.57X102  cm 

Conical  beam  width,  ©(• 

1.60  deg 

1.60  deg 

Antenna  gain,  G 

1.4X104 

8.0X103 

Minimum  detectable  signal 

1.6X10-14  w 

1.3X10" 14  w 

The  receiver  band-pass  characteristics  of  the  two  main  components 
have  been  matched  by  tuning  the  WSR/TAM-1  receiver  (originally  from 
an  AN/CPS-6  radar)  from  1.2  Mhz  to  1.0  Mhz  and  the  AN/CPS-9  receiver 
from  0.5  Mhz  to  1.0  Mhz,  both  on  long  pulse.  The  receivers  are 
calibrated  using  two  signal  generators:  TS-147D/UP  for  the  3.2-cm 
aet,  and  TS-155C/UP  for  the  10.3-cm  set.  Calibration  is  accomplished 
no  less  frequently  than  once  a  week,  more  frequently  as  operating 
conditions  permit. 
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An  unsuccessful  effort  was  made  to  acquire  from  the  National 
Science  Foundation  funds  to  procure  a  mono-pedeatal  capable  of 
supporting  a  single  array  of  three  antennas  with  the  centers  of  the 
dishes  at  the  same  height  above  the  ground.  It  now  appears  that  such 
acquisition  must  await  completion  of  the  new  11-story  Oceanography 
and  Meteorology  Building  that  has  been  approved  by  the  Board  of 
Directors  of  the  Texas  A&M  System.  The  pedestal  has  been  Included  as 
special  equipment  for  this  building,  which  it  is  hoped  will  be  ready 
for  occupancy  no  later  than  September  1971. 

Meteorological  conditions  near  the  antenna  site  are  being 
recorded  by  means  of  standard  instruments  mounted  on  the  AN/CPS-9 
tower  (wind  speed  and  direction)  and  on  the  observing  platform  near 
the  base  of  the  WSR/TAM-1  tower  (short-wave  radiation,  tR  radiation, 
temperature,  dew  point,  rainfall).  Hourly  airways  observations  are 
filed  from  the  Federal  Aviation  Agency  station  at  Easterwood  Field, 

2  miles  southvest  of  the  radar  site.  Conditions  representative  of 
the  Eist  Yegua  Creek  Watershed  will  be  determined  during  the  coming 
year  by  means  of  two  automatic  micrometeorologic'l  stations  recently 
installed  (see  Fig.  II1-1).  A  description  of  these  stations  is 
given  in  Section  IV  of  this  report. 

The  first  year  or.  the  work  has  been  devoted  largely  to  the 
collection  of  data.  By  chance,  the  eastern  portion  of  Texas  has 
been  subjected  to  unusually  heavy  amounts  of  precipitation  during 
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FIG.  m-  I  EAST  YEGUA  BASIN  NEAR  DIME  BOX,  TEX . 


this  period*,  and  consequently  project  personnel  have  been  kept  busy 
just  operating  the  radars.  Some  data  extraction  has  been  rccomplished, 
although  analysis  of  the  data  will  have  to  <*w*it  compilation  of 
meaningful  statistical  samples. 

The  Principal  Investigator  of  Sub-task  b  has  devoted  a  primary 
portion  of  his  research  time  to  an  investigation  of  an  observation  of 
unusual  "angel"  echoes  that  was  recorded  on  both  radars  on  May  23,  1968. 
(Examples  of  the  PPI  and  RHI  photographs  were  presented  in  Quarterly 
Progress  Report  t 3.)  These  observations  were  made  during  a  period  when 
an  instrumented  aircraft  from  the  National  Center  for  Atmospheric 
Research  was  flying  cloud  penetrations  for  the  Department  of  Meteorology. 
The  aircraft  was  sent  aloft  into  the  area  of  radar  returns;  a  straight 
climb-out  from  the  ground  was  conducted  to  10,000  feet,  with  a  deck  of 
8ttatocumulu6  being  encountered  between  2,800  and  5,500  feet.  No 
clouds  were  visible  above  this  layer.  A  spiral  ascent  then  was  made 
from  10,000  to  20,000  feet  in  the  area  of  radar  echo.  Numerous 
fragments  of  inversion  layers  (as  indicated  by  stratified  concentra¬ 
tions  of  haze)  were  visible  above  and  around  the  aircraft  on  all  sides. 
When  the  aircraft  reported  an  altimeter-indicated  altitude  of  18,000 
feet,  the  radar-indicated  elevation  was  12,000  feet  and  a  strongly 
enhanced  return  was  received  from  the  10-crn  transponder  on  the  aircraft. 

_ 

For  example,  the  climatological  rainfall  average  for  College  Station 
is  38.7  inches  per  annum;  by  July  9,  1968,  more  than  42  inches  had 
been  recorded  in  1968  —  with  almost  6  months  of  the  year  to  go  I 
This  is  Che  first  time  that  the  annual  average  had  been  exceeded 
before  the  end  of  July. 


55 


This  enhancement  of  signal  (aa  Indicated  on  the  A-acopa)  persisted 
throughout  the  flight. 

The  aircraft  was  outfitted  to  record  observations  of  static 
pressure,  temperature,  dew  point,  liquid-water  content,  true  air 
speed,  and  heading,  among  other  parameters.  Because  of  the  unique 
appearance  of  the  echoes  (no  euch  returns  previously  had  been  observed 
from  this  site),  the  Principal  Investigator  will  continue  to  devote 
his  attention  to  the  analysis  of  the  attendant  propagation  conditions. 
It  is  anticipated  that  the  assistance  of  a  graduate  student  will  be 
enlisted. 

* 

Drs.  Runnels,  Moyer,  and  Clark  collaborated  In  a  paper  that  was 
presented  at  the  13th  Radar  Meteorology  Conference,  Montreal,  Quebec, 
August  20-23,  1968.  Appropriate  acknowledgment  of  Project  Themis 
support  was  given. 

C.  Summary 

The  following  table  shows  some  operational  statistics  of  Sub¬ 
task  b  during  the  first  year. 

Table  III-2.  Operational  Statistics  1967-68 
Number  of  Storm  Days  Observed  136 

Number  of  Aircraft  Research  Flights  Tracked  21 

Hours  Spent  In  Collecting  Radar  Data  1,252 


* 

Runnels,  R.  C. ,  Moyer,  V.  £. ,  and  Clark,  R.  A.,  Uncertainties  in 
the  Measurement  of  Liquid-Water  Content  Using  the  Gain-Step  Technique," 
pp.  388-391,  Proceedings ,  13th  Weather  Radar  Conference,  Montreal, 
August  20-23,  1968. 
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Table  III-2— Continued 


Hour*  Spent  in  Tracking  Research  Aircraft  28 

Man-Houra  Spent  in  Support  of  Aircraft  Flights  364 

Man-Hours  Lpent  in  Calibrating  Radars  165 

Man  Hours  Spent  in  Maintaining  Radars  1,078 

Man-Hours  Spent  in  Maintaining  Rain-Gage  Network  218 

Man-Hours  Spent  in  Extracting  Radar  Data  241 

Man-Hours  Spent  in  Analysing  Radar  Data  554 

Feet  of  Radar  Film  Exposed  28 , 850 

Feet  of  Radar  Film  Developed  25,250 

Feet  of  Radar  Film  Edited  25,900 

Feet  of  Radar  Film  from  Which  Data  have  been  Extracted  615 

Feet  of  Radar  Film  Catalogued  49,800 

The  activities  for  the  second  year  will  necessarily  follow  a 
similar  pattern  with  Increasing  analyses  as  the  data  bank  approaches 
significant  levels. 
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‘IV.  MODIFICATION  AND  INSTALLATION  OF  AUTOMATIC  METEOROLOGICAL  STATIONS 

William  H.  Clayton 


A.  Station  Modification! 

Thaaa  station*,  designated  Station  A  and  Station  B,  vara  pre¬ 
viously  employed  In  Project  Green  Glow  and  the  Dallas  Tower  Network 
program  to  provide  measurement  of:  wind  speed,  wet  bulb  temperature, 
and  dry  bulb  temperature  at  1/4,  1/2,  1,  2,  4,  8,  16,  and  32  n;  wind 
direction  at  6  m;  soil  temperature  at  3,  6,  12,  25,  40,  65,  and  100  cm; 
insolation;  albedo;  net  radiation;  and  coll  heat  flux.  Recording  of 
these  data  is  according  to  the  format  shown  in  Table  IV-1,  The  term 
"automatic" ,  as  applied  to  these  stations,  designates  a  capability  for 
remote  operation,  plus  the  fact  that  all  sensor  outputs,  through 
suitable  function  generation  and  amplification,  are  recorded  in 
conventional  meteorological  parameter  form.  More  correctly,  the 
sensor  outputs  are  converted  to  DC  voltages  with  scale  factors  of  1, 

10,  or  100.  The  output,  in  addition  to  the  hard  copy  format  indicated 
in  Table  IV-1,  is  also  stored  on  punched  paper  tape. 

Inasmuch  as  the  design,  construction,  and  operation  of  these 
stations  has  been  fully  reviewed  ..n  another  report  (Clayton  and 
Eckelkamp,  1963),  discussion  here  will  be  limited  to  the  modifica¬ 
tions  incorporated  in  the  stations  for  Project  Themis  use  plus 
station  location  details. 

The  purpose  of  these  stations  in  the  current  research  program 
is  twofold.  The  first  of  these,  associated  with  Sub-task  b,  is  to 
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Table  "Y-l.  Printing  Format  for  Any  Five-Minute  Period 


peek  methods  to  correlate  radar  coverage  vlth  the  micrometeorologieal 
environment.  The  second,  associated  with  Sub-task  d,  is  to  provide 
measurements  on  a  mean  scale,  with  simultaneous  direct  measurements 
of  the  Reynold's  stresses,  from  other  facilities,  in  order  to 
investigate  postulated  concepts  concerning  equivalence  of  the  latter 
In  terms  of  the  former. 

The  modifications  effected  to  these  automatic  meteorological 
stations  prior  to  Installation  in  the  Yegua  Watershed  are  listed 
below. 

1.  Elimination  of  the  soil  heat  flux  measurements  in  favor 
of  surface  temperature  measurements. 

2.  Replacement  of  Beckman  Whitley  Radiometers  (Model  N188-01) 
by  Thomthwaite  Radiometers  (Model  605) . 

3.  Replacement  of  wind  speed  measurements  at  8  levels  through 
hard  tube  circuitry  and  stepping  switch  logic  with  solid  state 
components. 

Elimination  of  the  soil  heat  flux  measurements  is  primarily  a 
consequence  of  analyses  of  data  collected  at  Dallas,  Texas  and 
Richland,  Washington  that  shows  the  soil  heat  flux  is  rare.ly,  if 
ever,  accurately  measured  by  soil  heat  plates  (King,  lq63). 

The  soil  heat  flux  can  be  determined  from  measurements  of  soil  tem¬ 
perature  profiles,  which  are  already  incorporated  in  the  measurement 
complex;  however,  measurement  of  surface  temperature,  unless  made 
directly,  must  be  extrapolated  from  the  soil  and  air  profiles  which 
necessarily  degrades  the  flux  computations. 
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Another  requirement  also  entered  In  the  decision  to  replace 
soil  heat  flux  measurements  In  favor  of  measurements  of  surface 
temperature.  This  requirement  Is  based  on  the  following  equation 
employed  In  a  soil  model  under  study  in  Sub-task  c. 


q 


s,o 


(T  -  T')/G 
o  o 


1) 


This  equation,  first  introduced  by  Halstead,  et  al.  (1957), 
relates  the  soil  heat  flux  at  the  surface  to  the  surface  temperatures 
from  the  air  side  and  the  soil  side,  and  these  values  will  differ 
according  to  the  thermal  insulation  provided  by  vegetation  and  other 
surface  debris.  Virtually,  nothing  is  known  quantitatively  of  this 
factor  of  thermal  resistance  (G)  other  than  its  value  for  typical 
surfaces  which  appear  to  range  between  100  to  5000  cm  /sec  deg/cal. 
Analyses  of  collected  data,  from  these  stations  at  Dallas  and  Richland, 
to  evaluate  G  from  Equation  1)  have  been  something  less  than  useful. 
However,  assuming  that  overly  large  errors  were  not  made  in  these 
analyses  in  determining  T  and  T'  by  extrapolation,  and  q  by 

O  O  8  »© 

application  of  the  Fourier  conduction  equation;  it  would  appear  that 
G  is  not  a  constant  for  a  particular  surface  but  varies  rather  widely 
with  soil  moistness  and  air  temperature.  On  the  other  hand,  if  a 
mean  value  of  G  for  a  particular  surface  is  assumed  constant,  the 
scatter  of  Tq  and,  particularly,  values  are  within  reasonable 
limits  of  extrapolation  error.  Hopefully,  direct  measurement  of  Tq 
will  resolve  this  question. 
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Several  Mena  for  measuring  surface  temperature  were  investigated 
and  the  method  of  infrared  thermometry  (Barnes,  Model  IT-3)  was 
chosen.  The  sensor  is  mounted  on  the  tower  at  approximately  a  6  ft 
elevation  and  "sees"  the  surface  near  the  base  of  the  tower.  Inas¬ 
much  as  the  output  of  the  device  for  the  normal  temperature  ranges 
is  essentially  linear,  no  function  generation  is  employed  in 
converting  the  sensor  output  before  recording.  Also,  in  view  of 
the  relatively  high  cost  of  this  sensor  and  its,  as  yet,  improved 
capability  to  actually  measure  the  surface  temperature,  only  one 
unit  is  employed  (Station  A).  Subsequently,  assuming  operation  is 
satisfactory,  an  additional  unit  will  be  Installed  at  Station  R 
and  function  generation  will  be  Incorporated  in  the  measurement 
circuitry  to  compensate  for  the  nonlinearity  involved  between  -50 
and  1-50*  C.  Inasmuch  as  the  surface  temperature  measurement 
replaces  the  soil  heat  flux  measurement,  it  will  be  recorded  as 
Indicated  in  Table  IV-1;  that  is,  at  the  end  of  the  fourth  minute 
of  each  five  minute  Interval. 

Replacement  of  the  Beckman  Whitley  net  exchange  radiometer 
by  the  Thomthwalte  radiometer  is  based  on  testing  of  both  instru¬ 
ments  during  the  Dallas  Tower  Program.  In  the  tests  of  the  Beckman 
Whitley  device,  which  were  comparison  rather  than  calibration  tests 
inasmuch  as  the  latter  was  not  possible  to  project  facilities, 
five  radiometers  including  one  standard  (normally  maintained  in 
sealed  storage)  were  mounted  side  by  side  and  near  simultaneous 
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measurements  were  made  in  up-wind,  down-wind,  and  cross-wind 
configurations.  Also  measurements  in  an  Inverted  position  In  the 
same  wind  orientation  were  made.  These  tests  show  that  the  vari¬ 
ability  of  one  instrument  to  another  with  regards  to  tracking  capa¬ 
bilities  was  no  better  than  20%,  and  that  the  Instrument  response 
was  strongly  affected  by  wind  velocity.  Oddly  enough,  the  cross- 
wind  configuration  yielded  better  results  than  the  down-wind  config¬ 
uration  which  is  recommended  by  the  manufacturer.  The  Instrument 
response  was  not  significantly  varied  by  plate  reversal. 

Similar  comparison  with  two  Thornthwaite  radiometers  (supplied 
by  Dr.  L.  Frltschen  for  test  purposes)  showed  that  the  tracking 
capabilities  under  all  conditions  were  no  worse  than  5%  and,  of 
course,  the  sensors  were  unaffected  by  wind  flow. 

Modification  of  the  wind  measurements  system  applies  to  the 
conversion  of  the  sensor  output  rather  than  he  sensors  themselves. 
In  actuality  the  basic  correction  sought  in  this  system  was 
replacement  of  the  wind  counters  (Neuron  Type  7005D90AS)  by  more 
reliable  counters.  However,  investigation  of  several  other  counters 
did  not  demonstrate  significant  potential  for  improvement  in  this 
regard.  Since  counter  failure  often  occurred  in  the  units  digit 
(highest  operating  frequency),  replacement  of  the  five  digit  counter 
by  single  digit  counters  appeared  to  afford  operational  capability 
which  was  equivalent  to  an  improvement  in  counter  reliability,  in 
that  a  single  counter  failure  would  not  require  replacement  of  a 
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couplet*  five  digit  counter  as  had  been  the  caae  previously.  This 
modification,  because  of  the  radically  different  operating  character* 
Istics  of  the  single  counters,  introduced  the  necessity  for  several 
other  modifications  which  in  total  exceeded  ti.o  efforts  and  cost 
associated  with  the  design  of  a  completely  new  wind  measurement  system. 

In  the  following  description  of  the  modified  wind  system,  the 
discussion  is  patterned  after  that  given  in  the  previously  cited 
report  covering  the  operation  of  these  automatic  meteorological 
stations  and,  specifically,  replaces  pages  37-' 0  and  88-90; 

Figures  20,  21,  and  33;  and  Tables  4a,  4b,  and  18  in  that  report. 

The  hardware  modifications  to  the  wind  system  consists  of  nine 
types  of  printed  circuit  subassemblies  and  a  mainframe  that  holds 
these  assemblies  as  well  as  the  Veeder-Root  (Model  AG-196901)* 
decade  counters.  The  latter  directly  replaces  the  Neuron  counters 
of  the  old  wind  system.  Of  these  printed  circuit  assemblies,  seven 
are  removable  while  the  other  two  are  permanently  mounted  on  the 
front  panel.  Each  of  the  assemblies  is  identifiable  by  a  title 
corresponding  to  the  circuit  schematic  as  listed  in  the  following 
figures. 

Fig.  IV- 1.  7/11  and  Input  (8) 

Fig.  IV-2.  Reset  and  Set  20  Detector  (1) 

Fig.  IV-3.  Power  Driver  (2) 

Fig.  IV-4.  Lamp  Driver  (3) 

*The  counters  have  been  altered  internally  such  that  they  reset 
to  9  rather  than  9  or  0. 


Fig.  IV-5. 

Control  (1) 

Fig.  1V-6. 

Control  Relay  (1) 

Fig.  1V-7. 

Type  Relay  (1) 

Fig.  IV-8. 

Type  Control  I  (1) 

Fig.  IV-9. 

Type  Control  II  (1) 

The  numbers  in  parentheses  by  each  figure  title  indicates  the  number 
of  assemblies  employed  in  the  complete  system. 

The  signal  coding  employed  in  the  above  listed  figures  is 
described  in  Table  IV-2.  Use  of  this  coding  in  the  circuit 
schematic  descriptions,  which  follow,  would  likely  be  a  factor  of 
complication  and  is  avoided.  However,  the  correlation  to  the  signal 
coding  is  obvious 

a.  Fig.  IV-1  7/11  and  Input 

As  Implied  by  the  title,  this  circuit  receives  the  anemometer 
input  and  converts  the  anemometer  pulses  so  that  for  every  11  input 
pulses  there  are  seven  output  pulses.  This  is  accomplished  by  a 
scale-of-sixteen  binary  counting  circuit  that  resets  at  the  count 
of  eleven.  The  reset  is  accomplished  by  detecting  when  the  '1', 

'2',  and  '8'  outputs  are  all  logic  *1'  (5  volts).  The  three-input 
NAND  gate  which  detects  this  logic  state  then  sends  a  logic  'O' 
signal  (0  volts  or  absence  of  logic  '1')  to  another  three-input  NAND 
which  has  its  inputs  strapped  together.  The  second  NAND  then  supplies 
a  logic  '1*  to  a  third  and  final  three-input  NAND  whose  output  drives 
the  'direct  clear'  on  the  scale-of-sixteen  circuit.  The 
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Table  IV-2.  Signal  Coding 


SIGNAL  CODE 

LOCATION 

FUNCTION 

DIM 

TBI,  TB2 

DIGIT  1  READ  COMMAND 

D2R 

TB1,  TB2 

DIGIT  2  READ  COMMAND 

D3R 

TB1,  TB2 

DIGIT  3  READ  COMMAND 

D4R 

TB1 ,  TB2 

DIGIT  4  READ  COMMAND 

ED 

TB1,  TB2 

TYPING  ERROR  DETECT 

EDC 

TB1,  TB2 

ERROR  DETECTOR  RESET 

H 

TB1 ,  L7-19 

TYPE  HYPHEN  COMMAND 

HP 

TB2 ,  L7-4 

HIGH  PUNCH  COMMAND 

HPR 

TB1,  TB2 

HIGH  PUNCH  READ  COMMAND 

OSC 

TB2 ,  U10-2 

ERROR  DETECTOR  DRIVE 

RE 

TBI,  TB2 

RESET  ENABLE  COMMAND 
(AFTER  BUFFERS) 

REB 

TB1,  TB2 

RESET  ENABLE  COMMAND 
(BEFORE  BUFFERS) 

KS 

TBI,  U 10-14 

RESET  COMMAND 

RSF-1 

TBI,  Ul-18,  L5,-5 

RESET  FAIL  -  REGISTER  1 

RSF-2 

TBI,  U2-18,  L5-7 

RESET  FAIL  -  REGISTER  2 

RSF-3 

TBI,  U3-18,  L5-9 

RESET  FAIL  -  REGISTER  3 

RSF-4 

TBI,  U4-10,  L5-11 

RESET  FAIL  -  REGISTER  4 

RSF-5 

TB2,  U5-18,  L5-13 

RESET  FAIL  -  RECISTER  5 

RSF-6 

TB2 ,  U6-18,  L5-15 

RESET  FAIL  -  REGISTER  6 

RSF-7 


TB2 ,  U7-18,  L5-17 


RESET  FAIL  -  REGISTER  7 


Table  IV-2— - Continued 


SIGNAL  CODE  *  LOCATION  FUNCTION 


RSF-8  TB2,  U8-18,  L5-19  RESET  FAIL  -  REGISTER  8 


R1 

TB1 , 

TB2 

R2 

TB1, 

TB2 

R3 

TB1, 

TB2 

R4 

TB1, 

TB2 

R5 

TB1, 

TB2 

R6 

TB1, 

TB2 

R7 

TB1, 

TB2 

R8 

TB1, 

TB2 

R9 

TB1, 

TB2 

RO 

TB1, 

TB2 

R1R 

TB1 

R2R 

TB1 

R3R 

TB1 

R4R 

TB1 

R5R 

TB1, 

TB2 

R6R 

TB1, 

TB2 

R7R 

TB1 , 

TB2 

R8R 

TB1, 

TB2 

Rl-Dl 

TB1 

R1-D2 

TB1 

R1-D3 

TB1 

TYPE  '1'  ENABLE 

TYPE  '2'  ENABLE 

TYPE  '3'  ENABLE 

TYPE  '4'  ENABLE 

TYPE  '5'  ENABLE 

TYPE  '6'  ENABLE 

TYPE  '7'  ENABLE 

TYPE  '8'  ENABLE 

TYPE  '9'  ENABLE 

TYPE  'O'  ENABLE 

READ  REGISTER  1 

READ  REGISTER  2 

READ  REGISTER  3 

READ  REGISTER  4 

READ  REGISTER  5 

READ  REGISTER  6 

READ  RECISTER  7 

READ  REGISTER  8 

READ  REGISTER  1  -  DIGIT  1 

READ  REGISTER  1  -  DIGIT  2 

READ  REGISTER  1  -  DIGIT  3 
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1 

i 

i 

Tab  la 

IV-2 — Continued 

i 

i 

v 

j  SIGNAL  CODE 

LOCATION 

FUNCTION 

! 

R1-D4 

TB1 

READ  REGISTER  1 

-  DIGIT  4 

j  R2-D1 

I 

TB1 

READ  REGISTER  2 

-  DIGIT  1 

'  R2-D2 

1  [ 

TB1 

READ  REGISTER  2 

-  DIGIT  2 

1  R2-D3 

1 

TB1 

READ  REGISTER  2 

-  DIGIT  3  i 

|  . 

R2-D4 

TB1 

READ  REGISTER  2 

-  DIGIT  4  | 

R3-D1 

TB1 

READ  REGISTER  3 

-  DIGIT  1  | 

R3-D2 

TBl 

READ  REGISTER  3 

-  DIGIT  2  j 

R3-D3 

TB1 

READ  REGISTER  3 

-  DIGIT  3  ' 

R3-D4 

TBl 

READ  REGISTER  3 

-  DIGIT  4  ! 

i 

R4-D1 

TBl 

READ  REGISTER  4 

-  DIGIT  1 

!  R4-D2 

TBl 

READ  REGISTER  4 

-  DIGIT  2 

R4-D3 

TBl 

READ  REGISTER  4 

-  DIGIT  3  1  1 

i 

i  .  R4-D4 

TBl 

READ  REGISTER  4 

-  DIGIT  4  I 

1 

R5-D1 

TB2 

READ  REGISTER  5 

-  DIGIT  1  ' 

|  R5-D2 

TB2 

READ  REGISTER  5 

-  DIGIT  2  j  , 

* 

j  R5-D3 

TB2 

READ  REGISTER  5 

-  DIGIT  3 

j  R5-D4 

TB2 

READ  REGISTER  5 

-  DIGIT  4 

{ 

f  |  R6-D1 

TB2 

READ  REGISTER  6 

-  DIGIT  1 

R6-D2 

i 

TB2 

READ  REGISTER  6 

-  DIGIT  2 

R6-D3 

TB2 

READ  REGISTER  6 

-  DIGIT  3 

R6-D4 

*  ; 

! 

TB2 

READ  REGISTER  6 

-  DIGIT  4 

i 

i  • 

j 

• 

; 
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T«ble  IV-2 — Continued 


SIGNAL  CODE 

LOCATION 

FUNCTION 

R7-D1 

T62 

READ  REGISTER  7  -  DIGIT  1 

R7-D2 

TB2 

READ  REGISTER  7  -  DICIT  2 

R7-D3 

TB2 

READ  REGISTER  7  -  DIGIT  3 

R7-D4 

TB2 

READ  REGISTER  7  -  DIGIT  4 

R8-D1 

TB2 

READ  REGISTER  8  -  DIGIT  1 

R8-D2 

TB2 

READ  REGISTER  8  -  DIGIT  2 

R8-D3 

TB2 

READ  RECISTER  8  -  DIGIT  3 

R8-D4 

TB2 

READ  REGISTER  8  -  DIGIT  4 

S2F-I 

TB1,  Ul-17,  L5-5 

SET  20  FAIL  -  RECISTER  1 

S2F-2 

TB1,  U2-17 ,  L5-7 

SET  20  FAIL  -  REGISTER  2 

S2F-3 

TB1,  U3-17 ,  L5-9 

SET  20  FAIL  -  REGISTER  3 

S2F-4 

TB1,  U4-17 ,  L5-1I 

SET  20  FAIL  -  REGISTER  4 

S2F-5 

TB2 ,  U5-17,  L5-13 

SET  20  FAIL  -  REGISTER  5 

S2F-6 

TB2,  U6-17,  L5-15 

SET  20  FAIL  -  RECISTER  6 

S2F-7 

TB2,  U7-17,  L5-17 

SET  20  FAIL  -  REGISTER  7 

S2F-8 

TB2,  U8-17 ,  L5-19 

SET  20  FAIL  -  REGISTER  8 

T 

TB2,  L7-20 

TYPE  TABULATE  COMMAND 

TC 

TBI,  L7-18 

TYPE  COMPLETE  COMMAND 

TE 

TB1 ,  TB2 

TYPE  ENABLE 

TGR 

TBI,  TB2 

TYPE  TABULATE  ENABLE 

1 

TB2,  L6-20 

TYPE  ' 1 '  COMMAND 
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Table  IV- 2 — Continued 


SIGNAL  CODE 

LOCATION 

2 

TB2, 

L6-19 

TYPE 

3 

TB2, 

L6-18 

TYPE 

4 

TB2, 

L6-17 

TYPE 

5 

TB2 , 

L6-6 

TYPE 

6 

TB2 , 

L6-5 

TYPE 

7 

TB2 , 

L6-4 

TYPE 

8 

TB2, 

L6-3 

TYPE 

9 

TB2, 

L6-2 

TYPE 

0 

TB2, 

L6-21 

'  YPE 
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FUNCTION 
2'  COMMAND 
3*  COMMAND 
4'  COtMAfcD 
S'  COMMAND 
6'  COmAND 
7 1  COMMAND 
8'  COMMAND 
9'  COMMAND 
O'  COMMAND 


I 


V 


FROM 
TBI  -*XD  9 
FIG.  13 

MOM 

TBI  -*XD  0 
TBI  -*XD  2 
FIG. 13 


DETECT 
9  • 
INPUT 


I _ _  1/2  IC 


6  - 1  1/2  830P 


DETECT  0  \  5- 
INPUT  214- 


SET  20  OSC 


FROM  FIG.  5  — *  CONTROL  CARD 

SET  20  OSC  FAIL 

3 

FROM  FIG.  5—*  CONTROL  CARD 

RESET 

2 

FROM  FIG  5-*  CONTROL  CARD 

RESET  FAIL 

FROM  FIG.  5  — ►  CONTROL  CARD 


7/11  COUNTER 


COUNT  ON 

I 

FROM  FIG.  5— ►  CONTROL  CARD 


1/2  830  P 
1/2  IC  I 


R2  <2.2K 


jL  891  p  t 


cl  2.2  (TIMING  CAPACITOR) 


ANEMOMETER  INPUT 
5  FIN  HEX  PLUG 
FIG.  II 


^ — LOGIC  'o’  WHEN 
ANEMOMETER -FIRES 


40458  OR  4023) 


R,  <5.6  K 


x  •  i.E.s^.s.e.r.e 


FIG.  IV~I  INPUT 


NOTE:  (l)  NUMBERS  AT  THE  ENO  OF  A  SIGN/  LINE  OENOTE 
INPUT  OR  OUTPUT  PINS  C  ’  CIRCO.C  BOARD. 

(2)  NUMBERS  AROUND  COMPCUENTS  INDICATE 
CONNECTIONS  ON  IC'S. 


OUNTER 


.01 


« - 

R3 


IK  Rj 


_ _ _ . , _ ~ " 

_ >— !1L  Vli. 

ih-l|-  _yic 


IC  5 

♦Cp  859  P  Cf 


IC  6 

862  P 

3 

4l§/ 


»ING  CAR4CIT0R) 


l> 


{> 


IC  7 


836  P 


{> 

> 


RESET  OK 


TO  FIG  2 
-20 


RESET  AND  SET  20  DETECTOR 


RESET  ?0  OK  — •  TO  FIG.  2 
- 19 


RESET  AND  SET  20  OETECTOR 


12- 


SET  20  FAIL  LAMP  DRIVE 

- 17 

FIG.  12 

RESET  FAIL  LAMP  DRIVE 

- ie 

+  5  FIG.  12 


ADJUST! 

AiCf  TO  'S 


fp 


1C  8 
851  P 


ADJUSTS  PULSE  WIDTH 
Rg  TO  15  m  SEC. 

REGISTER 


ORIVE  -♦  TO 
-21 


Jh 


'SL 


POWER  DRIVER 
(SET  20  AND 
COUNT  DRIVE) 
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reason  for  this  cascading  is  to  allow  sufficient  delay  between  the 
fall  of  the  input  pulse,  which  caused  the  count  of  eleven  to  appear, 
and  the  clear  signal.  If  this  were  not  done  the  counter  might  clear 
and  enter  another  count  on  the  fall  of  the  same  pulse.  The  seven 
pulses  are  picked  off  through  two  R-C  networks  which  drive  two 
inverting  stages  whose  outputs  are  strapped  in  parallel.  This  com¬ 
bined  output  then  drives  a  monostable  circuit  (monostable)  that 
supplies  the  final  output  from  the  printed  circuit  assembly  (card). 
The  width  of  the  output  is  adjusted  to  15  msec  by  a  potentiometer 
located  on  the  card.  This  pot  is  available  from  the  front  of  the 
mainframe  cabinet  when  the  cabinet  is  open. 

From  the  diagram  it  is  apparent  that  the  input  pulse  drives 
the  scale-of-sixteen  circuit  through  a  four-input  gate  and  a  mono¬ 
stable.  The  purpose  of  the  monostable  is  to  shape  the  pulse  for 
reliable  counting.  The  gate  is  connected  to  three  lines  which  come 
from  the  'count  on'  line,  and  from  the  reset  outputs  of  the  dual 
flip-flop  circuit  (F/F)  located  on  the  card.  The  gate  is  enabled 
whenever  all  three  lines  have  logic  '1*  appearing  on  them.  When¬ 
ever  this  condition  is  met  the  output  of  the  gate  follows  the  input 
pulaes  from  the  anemometer.  The  'count  on'  line  is  energized  by  tbe 
control  card  and  the  F/F'e  are  in  the  reset  condition  whenever  the 
counters  reset  properly  after  the  preceding  RESET  signal.  The 
purpose  of  these  F/F's  is  to  detect  reset  and  'set  20'  errors  and 
in  the  event  of  no  errora  the  F/F's  remain  rsset.  In  the  event 
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that  failure  occurs  In  resetting,  either  one  or  both  F/F'a  will  set. 

If  this  happens,  an  error  lamp  will  light  on  the  front  panel  and  the 
typewriter  will  type  four  hyphens  to  indicate  failure.  Each  time 
a  RESET  is  received  from  the  external  main  system  the  wind  system 
circuitry  will  attempt  to  reset  all  counters.  Consequently,  counters 
which  are  intermittently  failing  to  reset  will  be  returned  to  the 
system  whenever  they  reset  properly  and  will  be  blocked  from  the 
system  whenever  they  fail  to  reset. 

Failure  to  reset  and  failure  to  'set  20'  signals  are  generated 
by  the  control  card  and  will  be  discussed  later. 

The  two-input  NAND  gate,  with  inputs  marked  'O'  Input  and  '2' 
input,  and  the  four-input  NAND,  with  inputs  marked  '9’s  input',  detect 
when  the  counter  assembly  associated  with  a  particular  input  card  has 
reset  to  a  reading  of  9999  and  when  it  has  set  a  reading  of  0020. 

The  outputs  of  these  gates  send  a  signal  through  an  inverting  stage 
to  the  system  reset  and  'set  20'  detector  card  and,  also,  either 
block  or  enable  the  reset  and  'set  20'  F/F  inputs. 

b.  Fig.  IV-2.  Reset  and  Set  20  Detector 
As  each  of  the  register  or  counter  assemblies  resets  to  a 
reading  of  9999  the  input  card  associated  with  a  particular  set  of 
registers  generates  a  logic  '1'  signal.  When  all  eight  register 
assemblies  have  completed  resetting, the  reset  detector  circuit 
generates  a  logic  '1'  signal  which  is  sent  to  the  control  card.  The 
control  card  then  changes  modes  of  operation  and  begins  to  set  0020 


on  the  register*.  As  each  register  arrives  at  0020  the  'set  20 
complete'  lines  from  each  input  card  go  to  logic  '1'  state.  When 
all  register  assemblies  read  0020  the  'set  20*  detector  sends  a 
logic  '1'  signal  to  the  control  card.  The  control  card  then  stops 
generating  the  'reset  20'  signal  and  generates  a  'reset  complete' 
signal  to  the  main  system. 

c.  Fig.  IV- 3.  Power  Drivers 

The  power  drivers  are  used  to  drive  the  count  and  reset  Inputs 
of  the  counter  assemblies.  In  the  case  of  the  count  inputs,  each 
of  the  eight  power  transistors  on  the  card  receives  pulses  from  the 
7/11  output  of  the  input  cards  and  actuates  the  count  input  of  the 
register  assesibly  associated  with  the  input  card  from  which  the 
pulses  were  received.  The  reset  inputs  receive  pulses  from  the 
second  power  driver  board  whose  input  comes  from  the  30  cycle  (5  volts) 
reset  oscillator  output  on  the  control  card. 

d.  Fig.  IV-A.  Lamp  Drivers 

The  lamp  driver  boards  receive  their  Inputs  from  three  sources: 
anemometer  input  'set  20  fall'  and  'reset  fall'.  The  lamp  driver 
in  the  anemometer  circuit  supplies  one  of  the  lamps  associated  with 
a  particular  register  assembly  to  indicate  whenevt.  vi  in;ut  pulse 
is  received.  The  'reset  fall’  lamp  driver  receives  s  inputs 
from  the  'reset  fail'  outputs  on  each  input  board  and  lights  a  lamp 
on  the  front  panel  in  case  of  RESET  failure.  The  'set  20  fell' 
lamp  driver  also  receives  its  Inputs  from  the  input  board  and  when 
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applicable,  lights  a  lamp  on  tha  front  panel  to  Indicate  that  a 
particular  reginter  assembly  has  failed  to  set  0020. 

e.  Fig.  IV-5.  Control 

With  the  exception  of  the  typing  format  the  control  card 
providaa  internal  lnattuctlons  to  the  wind  measurement  system.  Tha 
instructions  are;  'count  on',  reset,  'reset  oscillator',  'set  20 
oscillator',  'reset  fail',  and  'set  20  fail'.  A  30  cycle  square 
wave  is  also  provided  by  this  cird  as  well  as  a  50  volt  signal  back 
to  the  main  system  indicating  when  the  wind  counters  have  completed 
resetting.  The  control  card  receives  signals  from  both  the  counters 
and  the  main  system.  From  the  main  system  it  receives  50  volt 
pulses  which  instruct  tha  counter  assemblies  to  Btart  and  stop 
counting,  and  to  reset.  The  control  card  receives  logic  signals 
from  the  register  assemblies  and  'set  20*  detector  card  to  indicate 
whan  reset  and  'set  20'  instructions  have  been  completed.  Also,  a 
60  cycle  signal  is  fed  to  the  control  card  for  reduction  to  30  cycles. 
In  normal  operation  the  sequency  of  signals  would  be  start,  stop, 
reset,  'reset  complete'  and  'set  20  complete’.  On  receipt  of  a 
start  pulse  from  the  main  system  the  count  F/F  is  set.  The  count 
F/F,  when  set,  supplies  a  logic  '1'  signal  to  all  input  cards  to 
enable  the  input  gate.  Counting  in  all  registers  then  proceeds  until 
a  stop  pulse  is  received  from  the  external  mnin  system.  Since  the 
typing  cards  function  independently  they  will  be  discussed  later. 
After  typing  la  complete  the  main  system  generates  a  50  volt  RESET 
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command  which  resets  the  failure  F/F's  on  the  Input  boards  and 
resets  the  typing  boards.  When  the  reset  F/F  is  set,  a  two- Input 
NAND  gate  Is  enabled  which  allows  the  30  cycle  signal  to  be  applied 
to  an  inverter  which  drives  a  scale-of-alxteen  counter  and  the 
Inputs  of  the  reset  power  drivers.  When  all  registers  are  reset, 
the  'reset  complete*  signal  resets  the  F/F  to  remove  the  30  cycle 
output.  The  function  of  the  scale-of-slxteen  counter  Is  to 
determine  when  12  pulses  (30  cycle)  have  been  supplied  to  the  reset 
drivers.  If  12  pulses  are  supplied  without  a  'reset  complete' 
signal,  the  counter  Interprets  this  as  an  indication  of  reset 
failure,  since  a  maximum  of  only  nine  pulses  are  required  to  reset 
the  system,  and  generates  a  'reset  fail*  signal  of  logic  '1'  through 
the  two-input  HAND  gate  and  inverter  and  holds  this  signal  until 
a  count  (30  cycle)  of  IS  Is  received.  At  a  count  of  15  the  counter 
sends  a  logic  'O'  signal  back  from  the  four-input  NAND  gate  to 
reset  the  reset  F/F  and  clear  the  counter.  When  this  occurs  the 
'reset  oscillator'  and  'reset  fall'  signals  are  removed.  Further, 
through  the  combination  of  a  R-C  network  and  Inverter  connected  to 
the  reset  output  of  the  reset  F/F,  a  signal  is  generated  to  set  the 
'set  20'  F/F.  This  In  turn,  enables  a  two-input  NAND  to  permit 
transmission  of  the  30  cycle  signal  to  the  input  cards  and  to  a  F/F. 
The  Input  cards  then  proceed  to  count  21  pulses  (48  volts)  Into  the 
register  assemblies  (the  number  required  to  go  from  9999  to  0020). 

As  each  register  sets  0020,  the  Input  card  prevents  further  pulses 
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from  being  sent  to  the  register  by  disabling  a  two- input  HAND  gate 
on  the  input  board.  When  all  registers  have  set  0020,  a  logic  '1' 
signal  is  sent  from  the  reset  and  'set  20'  detector  card  which 
resets  the  'set  20'  F/F.  The  output  of  the  F/F  which  was  being 
driven  by  the  'set  20'  oscillator  signal  drives  the  input  of  a 
second  scale-of-sixteen  counter.  This  counter  is  connected  in 
In  the  same  way  as  the  one  associated  with  the  reset  circuit 
excepe  that  it  counts  only  every  other  pulse  because  of  the  F/F 
preceding  it.  This  allows  24  'set  20*  pulses  to  occur  before 
a  failure  signal  is  generated.  A  set  20'  failure  is  held  for  six 
counts  rather  chan  for  three  as  in  the  case  of  the  reset  F/F.  When 
the  'set  20’  F/F  is  reset  it  generates  a  logic  '1*  output  through 
the  R-C  and  inverter  combination  which  indicates  to  the  main  system 
that  the  total  wind  system  has  been  reset. 

f.  Fig,  IV-6.  Control  Relays 

The  control  relays  are  required  to  transform  the  logic  levels 
(5  volts)  used  in  the  wind  measurement  system  to  the  50  volt  level 
used  by  the  main  system.  There  are  six  relays  and  driving  circuits 
on  this  card  and  they  supply  signals  for  the  following  functions: 

1.  system  reset  complete 

2 .  typewriter  hyphen 

3.  typewriter  tabulate 

4.  high  punch  on  tape 

5.  type  wind  complete 

6.  tape  leader  advance 
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The  input*  for  tech  relay  driving  circuit  are: 

1.  system  reset  complete  from  control  card 

2.  typewriter  hyphen  from  type  control  I 

3.  typewriter  tabulate  from  type  control  II 

4.  high  punch  from  type  control  II 

5.  type  wind  complete  from  type  control  I 

6.  tape  leader  advance  from  reset  oscillator  output 

g.  Fig.  IV- 7.  Type  Relays 

The  type  relays  are  used  to  type  and  punch  the  digits  0  through 

9  and  receive  their  inputs  from  type  control  II. 

h.  Fig.  IV-8.  and  IV-9.  Type  Control  I  and  II 

The  type  control  system  contains  two  shift  registers  which 
are  used  in  conjunction  to  identify  the  particular  digit  to  be 
typed.  One  of  the  shift  registers  is  a  five-bit  unit  which  is 
closed  on  itself  while  the  other  register  is  an  eight-bit  unit 
which  is  open  at  the  end.  The  five-bit  register  is  used  to  identify 
the  digit  to  be  read  while  the  eight-bit  register  is  used  to 
identify  the  counter  assembly.  These  registers  make  up  the  bssic 
logic  of  the  system  while  the  balance  of  the  logic  Is  used  to  decode 
the  outputs  of  the  shift  registers  and  to  supply  the  necessary 
control  functions. 

On  receipt  of  a  ’type  wind*  command,  the  type  F/F  is  set  and 
the  two  shift  registers  have  a  logic  '1'  set  into  the  first  bit 
location.  When  the  type  F/F  is  set,  the  type  AND  gate  is  enabled 
and  a  logic  *  1'  is  applied  to  the  digit  AND  gates.  With  the  logic 
’1’  In  the  first  bit  locations  of  the  shift  registers,  a  logic  '1' 
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signal  will  k  ifflM  to  tki  op— oa  11m  of  the  tint  digit  U 
the  first  vlad  counter  assembly.  Depending  on  the  reeding  comtataad 
In  this  assembly,  this  logic  '1*  eigne  1  will  appear  o»  one  of  the 
digit  ARC  gate*.  This  signal,  together  with  the  signal  free  the 
type  AMD  gate,  cause*  one  of  the  'type  digit'  relays  te  be  e—rglaed. 
When  a  'typa  digit  relay  la  energised,  a  50  volt  elgnal  la  applied 
to  an  ootpet  11m  to  drive  the  appropriate  solenoid  on  the  type¬ 
writer.  On  completion  of  typing  a  digit,  the  typewriter  return*  a 
50  volt  typewriter  Interlock  signal  to  the  typo  control  beard,  title 
signal  causae  the  monos table,  which  Is  connected  to  the  type  AMD 
gate,  to  drive  to  a  logic  *0'  scats,  this  tonovea  the  logic  *1' 
signal  fron  the  digit  AMD  gate*  and  shifts  the  logic  '1*  signal  la 
the  first  bit  location  of  the  five-bit  register  to  th*  second  bit 
location.  After  the  aonoe table  has  dal a/ad  for  approximately 
35  Mac*  It  returns  to  logic  .  The  digit  AMD  gate*  arc  agalr 
enabled,  and  the  second  digit  la  the  first  counter  assembly  le 
typed.  After  all  four  digits  are  typed,  the  fifth  bit  la  the  fire- 
bit  register  energises  the  tab  AND  gats  and  the  typewriter  tabulate* 
to  the  nest  position.  On  receipt  of  the  typewriter  interlock  signal 
aftsr  tabulation  tbs  five-bit  register  again  has  a  logic  '1*  in  the 
first  bit  location.  When  the  fifth  bit  goes  from  logic  '1'  to  logic 
'O'  It  advances  the  logic  '1*  In  ehe  first  bit  location  of  the  eight- 
bit  register  to  the  second  bit  location.  This  sequencing  of  digits 
and  registers  continues  vntil  the  fourth  digit  in  the  eight 
counter  raasmbly  has  been  typed.  Ac  this  tine,  the  tab  signal 
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Is  disabled  r.ud  a  type  complete  signal  is  generated  which  sends  a 
SO  volt  pulse  back  to  the  main  system. 

In  the  event  a  digit  should  fail  to  type  after  f  \e  digit  AND 
gates  have  been  energized,  Che  system  generates  a  command  to  type 
a  hyphen.  This  command  Is  generated  by  a  scale-of-sixteen  circuit 
which  begins  counting  a  30  cycle  signal  sr  soon  as  the  digit  AND 
gates  are  enabled.  If  the  count  reaches  15  without  a  typewriter 
Interlock  signal  being  received,  the  scale-of-sixteen  circuit 
generates  the  type  hyphen  comand.  As  each  digit  is  typed  the 
scale-of-sixteen  circuit  is  reset  so  that  each  digit  is  allowed 
a  half  second  in  which  to  type  before  a  hyphen  is  typed  in  its 
place.  This  failure  circuit  is  also  used  whenever  a  counter  has 
failed  to  reset  properly.  If  reset  or  "set  20'  failure  occurs,  the 
output  from  the  eight-bit  register,  associated  with  the  particular 
counter  assembly,  is  disabled  so  that  when  the  data  in  that 
register  is  to  be  typed,  the  common  lines  to  each  of  its  digits 
are  not  energized  and  therefore  cannot  be  typed.  As  a  result  this 
is  Interpreted  as  four  successive  failures  and  the  system  types 
four  hyphens.  Thus,  a  typed  group  of  four  hyphens  indicates  that 
a  particular  counter  assembly  has  failed  to  reset  properly  or  has 
failed  (all  digits)  to  read  out  properly. 

All  of  the  preceding  references  to  typewriter  control  apply 
also  to  simultaneous  punched  paper  tape  commands. 
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Table  IV-3.  Wind  Speed  Cabling 


Connector 

Pin  # 

Function 

Connects  to 

A 

Photocell  supply 
voltage  +5  VDC 

Anemometer,  Pin  2 

A-l 

B 

Photocell  input 

Anemometer,  Pin  3 

through 

D 

3  VAC 

Anemometer,  Pin  4 

A- 8 

E 

Ground 

Anemometer,  Pin  1 

H 

Ground 

Anemometer,  Pin  1 

1 

Type  Wind 

Main  programmer.  Pin  W 

2 

Wind  typed 

Main  programmer ,  Pin  N 

3 

Beset 

Main  programmer.  Pin  a; 

timer  Pin  J 

C-l 

4 

Counter  reset 

Main  programmer.  Pin  Z 

5 

Start  counters 

Timer,  Pin  B 

6 

NC 

7 

Stop  counters 

Timer,  Pin  C 

8 

NC 

1 

115  VAC 

2 

+50  VDC 

3 

T.  I. 

C-2 

4 

Tape  advance 

Typer/punch  Junction  box. 

10  Pin,  Pin  I 

5 

50  VDC  return 

6 

115  VAC 

Table  IV-3. — Continued 


Connector 


C-3 


Pin 

t  Function 

Connects  to 

1 

Typewriter  punch 
junction  box 

15-Pin,  Pin  1 

2 

Typewriter  punch 
junction  box 

Pin  2 

3 

Typewriter  punch 
junction  box 

Pin  3 

4 

Typewriter  punch 
junction  box 

Pin  4 

5 

Typewriter  punch 
junction  box 

Pin  5 

6 

Typewriter  punch 
junction  box 

Pin  6 

7 

Typewriter  punch 
junction  box 

Pin  7 

8 

Typewriter  punch 
junction  box 

Pin  8 

9 

Typewriter  punch 
junction  box 

Pin  9 

10 

Typewriter  punch 
junction  box 

Pin  10 

11 

Typewriter  punch 
junction  box 

Pin  11 

12 

Typewriter  punch 
junction  box 

Pin  12 

13 

Typewriter  punch 
junction  box 

Pin  13 

14 

Typewriter  punch 
junction  box 

Pin  14 

15 

Typewriter  punch 
junction  box 

Pin  15 
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2  INPUT  OR  ( I809P) 


2  INPUT  NAND  (846P) 


3  INPUT  .NAND  (862P) 


4  INPUT  NAND  (830P) 


2  INPUT  NOR  ( I8I0P) 


^>o -  INVERTER  (836P) 


flip-flop,  monostable,  or  counter 

FLIP-FLOP’  852 P,  853P,  OR  848P 
MONOSTABLE'  851 P 
COUNTER’  839P 
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Figures  IV-10  through  IV-14  show  the  signal  flow  pattern  In 
the  modified  wind  system  and  circuit  wiring.  It  should  be  noted 
that  Integrated  circuit  (IC)  modules  in  these  and  previous 
schematic  diagrams  are  manufactured  by  Motorola.  The  IC  signal 
levels  are  5  volts. 

The  interconnection  of  the  modified  wind  system  to  the  control 
systems  employed  in  the  station  trailers  is  given  in  Table  IV-3. 

5.  Station  Installation  Characteristics 

The  station  locations  (see  Fig.  I1I-1)  and  configurations  are 
aeh«r at leal ly  represented  in  Figures  IV-15  and  IV-16.  The  areal 
extent  In  the  letter  figures  is  a  circle  of  500  ft  radius  centered 
at  the  InatnsMnt  tower. 

The  following  tables  show  the  soil  characteristics  at  these 
stselona. 


Table  IV-4.  Soil  Characteristics  st  Station  A  (Burleson  County  - 


30*  29  *1*  N  - 

96*  53  min  W) 

Depth 

1* 

0.  M. 
X 

t 

lb/acre 

lb /acre 

Ca 

lb/acra 

Soluble 

salts 

lb/acre 

Magnesium 

lb/acre 

Na 

lb/acre 

3  c* 

b.l 

1.3 

1.2 

116 

1060 

120 

140 

112 

6  cm 

6.3 

e.8 

2.4 

46 

620 

120 

110 

112 

12  cm 

6.1 

0.6 

1.2 

SO 

600 

i:o 

105 

100 

cm 

e.3 

C.5 

3.6 

30 

420 

120 

110 

100 

<*0  c« 

6.6 

0.6 

0 

13< 

1200 

120 

475 

100 

65  cm 

6.5 

0.5 

0 

304 

2540 

120 

>500 

240 

160  cs 

6.1 

0.2 

0 

324 

1220 

120 

>500 

476 

96 
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Table  IV-S.  Soli  Character is ties  at  Station  B  (Lee  County  30°  32  min  N  - 
96*  58  min  W) 


Depth 

PH 

0.  M. 
X 

P 

lb /acre 

K 

lb/acre 

Ca 

lb/acre 

Soluble 

salts 

lb/acre 

Magnesium 
lb /acre 

Na 

lb/acre 

3  cm 

6.7 

0.5 

6.2 

382 

2300 

384 

>500 

620 

6  cm 

6.9 

0.6 

10.1 

98 

940 

384 

170 

312 

12  cm 

7.7 

0.4 

20 

464 

2460 

1104 

265 

430 

25  cm 

5.2 

0.5 

1.2 

286 

2840 

288 

>500 

620 

40  cm 

5.4 

0.5 

1.2 

232 

2760 

432 

>500 

1134 

65  cm 

6.1 

0.4 

1.2 

196 

3800 

624 

>500 

1310 

100  cm 

7.1 

0.3 

1.2 

186 

3740 

624 

>500 

1630 

The  soil  type  at  Station  A  is  similar  to  that  classified  as  Bonham 
Loam  while  that  at  Station  B  Is  similar  to  the  classification  Tarrant 
Stoney  Clay.  Conductivity  and  specific  heat  determinations  have 
not  been  made  at  this  writing. 

Both  stations  are  located  in  pasture  land  In  sparsely  settled 
regions.  The  terrain  can  best  be  described  y  rolling  with 

areas  of  medium  to  dense  forestation  ranging  from  20  to  60  ft  in 
height.  These  forested  areas  are  scattered  and  tend  to  decrease  in 
extent  and  frequency  with  distance  from  the  Yegua  Creek. 

Stations  A  and  B  are  completely  self  contained  with  respect  to 
sensors,  measuring  systems,  maintenance,  calibration,  and  supporting 
equipment.  Power  is  supplied  from  commercial  sources  (R.  E.  A.)  as 
230-115  volts,  60  cycle,  Bingle  phase.  Operation  can  be  from  motor 
driven  generators. 
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As  previously  noted,  the  sensors  outputs  are  converted  to  DC 
voltages  and  so  scaled  as  to  be  numerically  equivalent  to  conventional 
meteorological  parameters.  Calibration  procedures  for  the  sensors 
and  the  digital  conversion  systea  can  be  routinely  performed  without 
modification  of  the  system  itself,  and  in-line  check  procedures  are 
Incorporated  to  permit  distinction  between  signal  indication  and 
correct  operation. 

The  trailers  housing  the  stations  are  10  ton  military  semi¬ 
trailer  vans  (U.  S.  Army  Model  SKD  3809)  of  approximate  dimensions 
26  x  8  x  6.2  ft  with  an  overall  height  of  approximately  13  ft.  The 
vans  have  been  modified  by  installation  of  five  windows,  a  side  door, 
an  air  conditioner  port,  an  interior  lining  of  1/2  in  plywood,  a 
floor  covering  of  asphalt  tile,  and  carrying  boxes  welded  to  the 
underside  of  the  trailer  for  transport  of  cables,  sensors,  and 
station  supplies. 

The  measurement  cower  is  of  triangular  steel  construction, 

16  in  on  a  side,  made  up  of  eleven  10  ft  sections  and  one  5  ft 
section  for  an  overall  height  of  115  ft  (Rohn  Model  40C).  Trian¬ 
gular  guying  at  30  ft,  65  ft,  and  95  ft  is  used  to  support  the 
tower  with  concrete  'dead  men'  located  75  ft  from  the  tower  base 
as  the  guy  anchors.  The  guys  themselves  are  3/8  in  aircraft  cable. 
The  tower  is  supported  by  a  concrete  block,  approximately 
2x2x4  ft,  in  which  the  necessary  bolts  for  securing  the  tower 
base  platform  are  embedded.  Erection  of  the  tower  is  accomplished 
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by  *  floating  gin  pole  which  damps  to  successive  tower  sections  as 
construction  proceeds.  The  weight  of  the  10  ft  tower  section  is 
approximately  65  lb.  The  instrument  arms  are  made  of  aluminum  channel 
and  3/4  in  aluminum  tubing  with  the  channel  bolted  to  the  tower  side 
aa  support  for  the  tubing  which  can  be  moved  in  and  out  from  the 
tower.  The  arm  installation  at  any  level  is  so  arranged  that  the 
anemometer  cups  and  thermocouples  are  at  the  designated  sensing 
elevation  with  the  anemometer  attached  to  the  free  end  of  the 
tubing,  5  ft  from  the  tower  side,  and  the  thermochelter  attached 
to  the  channel,  18  in  from  the  tower  side. 

This  completes  the  discussion  of  the  modification  and 
installation  of  the  automatic  meteorological  station  in  the  Yegua 
Watershed,  except  to  note  that  the  operation  of  these  stations  is 
conducted  under  Sub-task  b.  Operational  scheduling  and  analytical 
procedures  will  be  outlined  in  a  technical  report  of  that  sub- 

i  1 


task  at  a  later  date. 


REFERENCES 


Clayton,  W.  B. ,  and  B.  J.  Eckalkaap,  1963:  Inprovonent,  avaluatlon, 
and  ex tana ion  of  capabllltlea  of  meteorological  a Inula t ion. 
Final  Rot.  Contract  DA36-039  SC84942.  Texax  A  &  M  Reaaarch 
Fdn. ,  63-1T. 

Halstead,  N.  H. ,  R.  L.  Rlchnan,  W.  G.  Covey,  and  J.  D.  Merryman,  1957 
A  preliminary  report  on  the  deBign  of  a  computer  for 
micrometeorology.  J.  of  Met. .  14. 

King,  B.  L.  D. ,  1963:  An  evaluation  of  heat  flow  tranaducera  as 
a  means  of  determining  soil  heat  flow.  N.  Sc.  thesis, 

Texas  A&M  University. 


V.  ANALOG  MODELING 

William  H.  Clayton  and  F.  L.  Deacon 


The  Themis  portion  of  the  analog  modeling  research  being 
conducted  at  Texas  A&M  is  primarily  concerned  with  investigation 
of  the  applicability  of  various  meteorological  relationships 
commonly  employed  in  model  configurations.  During  the  current 
reporting  period,  the  main  emphasis  has  been  devoted  to  the 
frictional  boundary  layer  and  particularly  to  the  concept  and 
employment  of  exchange  coef ticients  to  provide  necessary  coupling 
of  atmospheric  model  layers. 

Although  exchange  coefficients  are  widely  accepted  bb 
physically  sound  parameters  describing  turbulent  flow  character¬ 
istics,  very  little  quantitative  information  is  available  and 
functional  definition  exists  only  as  hypotheses.  The  basic 
definition  of  the  exchange  coefficient  rests  in  the  following 
equation,  postulated  by  Wilhelm  Schmidt  (Sutton,  1953),  based  on 
an  analogy  of  turbulent  and  laminar  flow  conditions  as  applicable 
to  the  flux  of  momentum. 


dS 

pK  T 
m  dz 


1) 


The  term  K  in  this  relationship  is  the  eddy  coefficient 
tn 

for  momentum,  or  eddy  viscosity,  which  is  presumed  to  be  analagous 

to  the  kinematic  viscosity  of  a  fluid  in  laminar  motion.  The 

product  pK  is  the  exchange  coefficient, 
m 

The  concept  of  analogy  between  the  two  flow  patterns  may 
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or  nay  not  be  realistic.  Kr>'.  aver,  if  one  does  accept  the  analogy 
then  the  definition  of  the  exchange  coefficient  for  momentum  la 
as  defined  in  Equation  1) .  Similar  analagous  reasoning  leads 
to  corresponding  definitions  for  the  eddy  coefficients  for  heat 
and  vapor,  and  no  difficulties  are  encountered  until  an  attempt 
is  made  to  employ  these  coefficients  in  meaningful  fashion  in 
meteorological  models.  Such  application  necessarily  involves 
usage  of  other  meteorological  relationships,  often  of  the  widely 
accepted  variety,  and  this  can  lead  to  results  that  are  question¬ 
able  at  best. 

An  illustration  of  this  is  afforded  through  any  of  the  various 
diabatic  profiles  of  wind,  temperature,  and  vapor  pressure  that 
are  relatively  well  known  and  employed.  As  an  example,  consider 
the  following  relations  for  wind  and  temperature  given  by  Monin 
and  Obukhov  (1954) , 


dS 

dz 


kl  Vl> 


de 

oz 


kz^2^L 


2) 

3) 


where 

u*^0pC 

L  -  _  — - E. 

kgH 

and 

6  - - a— 

*  pc  ku* 

P  * 


4) 

5) 


All  of  the  symbols  employed  have  the  usual  meteorological 
designations  as  listed  in  the  glossary  at  the  end  of  this 


104 


report  section.  The  functions  and  are  not  uniquely  specified, 
but  if  the  eddy  coefficients  for  n omen tun  and  heat  are  assumed 
equal,  then  equal  Purther,  if  |s/l|<1,  and  s  »  *o,  then 
the  function  can  be  approximated  by 

*(f)  -  1  +  af  6) 

ignoring  terms  of  order  two  and  higher  in  a  series  expansion. 

Substitution  of  6)  in  2)  and  3)  leads  to  the  well  known 
Monin-Obukhov  profile  relationships. 

7> 

o 


The  value  of  a  suggested  by  Monin  and  Obukhov  is  .6.  Panofsky, 
Blackadar,  and  McVehll  (I960),  suggest  a  to  be  4.5.  No  matter 
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the  value,  other  than  the  trivial  aero,  it  is  seen  that  10)  can 

lead  to  negative  values  of  K  .  Also,  it  can  lead  to  K  values 

m  m 

that  Increase  with  decreasing  wind  speed.  Neither  of  these 
possibilities  are  physically  reasonable.  To  illustrate  the 
preceding,  Tables  V-l  and  V-2  below  show  K  .  values  obtained 

&i|h 

from  10)  for  various  wind  and  temperature  values  when  h  ■  8  m, 

—  2 

k  »  .4,  6  ■  300  #K,  g  ■  980  cm/sec  ,  and  a  ■  2  cm.  In  Table  V-l, 

o 

a  is  chosen  as  .6  whereat  in  Table  V-2,  a  is  chosen  as  4.5. 

Table  V-l.  K  Q  ValueB  (cra2/sec) 

0k|O 

a  ■  .6 


Sgin/see 

<es- 

e0rc 

8 

4 

0 

-4 

-8 

10 

19152 

20210 

21364 

22625 

24009 

5 

7134 

6638 

10682 

13612 

18147 

2.5 

1420 

2540 

5341 

21672 

-19096 

Table 

V"2*  Km,8 

Values 

(cm2/sec) 

a  ■  4.5 

SgOn/sec) 

<68- 

e0)*c 

8 

4 

0 

-4 

-8 

10 

10576 

14600 

21364 

34898 

74721 

5 

1144 

3040 

10682 

-47795 

-12406 

2.5 

118 

52 

5341 

-  3849 

-  3520 
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This  same  physical  inconsistency  is  met  whatever  pro tils 
relationships  are  usad.  However,  use  of  a  power  profile  for  the 
wl id,  and  omission  of  direct  consideration  of  a  dlabatic  temp¬ 
erature  profile  improves  the  picture  somewhat.  For  example,  use 
of  the  Deacon  profile. 


dS 

dr 


U*E 


-e 


kz 


(1  -  fO 


11) 


leads  to  the  following  expression  foT  the  eddy  viscosity  at 
height  z 


m 


Sk2(l  -  6)zo{1  "  6)zB 

t<f->a  “  B)  -  1] 


12) 


Use  of  12)  depends  on  a  functional  definition  of  0.  Professor 
Deacon  suggests  the  following  empirical  relationship 


0  -  1  -  1.43R1  -  6Ri2  -  10R13 


13) 


where  R1  is  the  layer  Richardson  number  between  the  surface  and 

elevation  h(z  ■  2  cm), 
o 

Substituting  13)  in  12)  and  using  the  same  values  of  h, 

0,  g,  and  k  employed  previously,  leads  to  the  K  values  shown  in 

in 

Table  V-3.  As  noted  previously,  all  inconsistency  is  not  removed  as 
can  be  seen  for  Sg  •  5  m/sec  and  ”  ®0  "  “®#c* 


107 


Sg(i»/aec) 


Tabic  V-3.  K  a  Values  (an  /sec) 

ra,o 


<e8 

~  v <C 

8 

4 

0 

-4 

-8 

15825 

18513 

21364 

24337 

27392 

2197 

5528 

10682 

16787 

22559 

0 

60 

5341 

15218 

26950 

In  the  discussion  of  the  Monin  and  Obukhov  relationship,  the 
question  of  physical  contradiction  rests  with  two  concepts;  the  log- 
linear  profiles,  and  the  eddy  viscosity  definition.  The  Deacon 
profile  introduces  another  factor,  Che  Richardson  number. 

In  view  of  the  extensive  employment  of  the  latter  as  a 
stability  criterion,  it  would  seem  illogical  to  question  its  worth. 
However,  a  careful  analysis  of  the  concepts  utilized  in  its  deri¬ 
vation  provides  some  understanding  why  its  determinations  from 
profile  measurements  leads  to  numerous  inconsistencies.  In  fact, 
one  is  led  to  believe  that  extensive  employment  of  the  Richardson 
number  could  be  more  a  matter  of  discussion  than  practice. 

Use  of  the  Richardson  number  in  objective  modeling  experiments 
requires  application  to  finite  layers.  Thu9  in  the  determination 
of  Table  V-3,  the  layer  Richardson  number  was  defined  as 

H  -  gia<V  V  141 

°  “S8 

Obviously  the  application  of  14)  to  light  wind  gradients  is 
limited.  Presumably,  one  can  retreat  from  such  situations  through 
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speculation  concerning  data  reliability  under  such  conditions,  but 
It  seems  feasible  to  question  also  the  Inherent  validity  of  a 
relation  that  is  applicable  to  continuous  phenomena  when  the  relation 
incorporates  discontinuity  (as  a  second  order  pole) . 

Several  possibilities  exist  to  rectify  this  difficulty  and 
the  addition  of  a  threshold  wind,  or  turbulent  anomaly  to  the  mean 
wind,  is  efficacious.  With  this  approach  14)  becomes 


Sz8<68  "  V 

RI  ■  9 

e(sg  +  b r 


15) 


Use  of  this  relation  (b  *  3  m/sec)  with  12)  and  13)  gives  the 

K  values  shown  in  Table  V-4. 
in 


Sg (m/sec) 

Table 

V'4-  Km,8 
(68- 

2 

Values  (cm  /sec) 

vc 

8 

4 

0 

-4 

-8 

20 

40540 

41630 

42727 

43835 

44951 

10 

18007 

19660 

21364 

23111 

24892 

5 

6517 

8489 

10682 

13024 

15436 

2.5 

1565 

3153 

5341 

7858 

10332 

1.25 

217 

1004 

?<*’9 

4777 

6547 

.63 

21 

304 

1 

2755 

3733 

.31 

2 

100 

668 

1498 

1998 

The  range 

of  modified  fti  in  Table 

V-4  is  0 

to  1.2. 
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None  of  the  values  given  in  Table  V-4  are  mathematically  or 


meteorologically  inconsistent  and  the  values  of  K  a~e  comparable 
to  those  obtained  with  other  profile  relations.  Further,  model 
solutions  of  the  parameters  of  the  frictional  boundary  layer  are 
mathematically  possible  for  wind  and  temperature  values  that  were 
previously  excluded  and,  no  model  degradation  of  previously 
included  values  occurs.  Of  course,  this  cannot  be  considered  as  an 
indication  of  meteorological  correctness  and  further  investigation 
is  underway.  Rather,  it  permits  model  experiments  for  further 
improvement  of  objective  f  recasting  in  the  frictional  boundary 
layer. 

No  rigorous  justification  exists  for  the  introduction  of  the 

mean  wind  speed  anomaly  in  the  layer  Richardson  number,  as  shown 

in  15),  other  than  its  usage  removes  mathematical  discontinuity. 

This  type  of  Justification  is  not  without  precedence  in  meteorology 

as  is  evidenced  by  the  term  z  in  tne  adiabatic  and  diabatic  profile 

o 

equations.  However,  zq  has  been  used  so  often  in  the  literature 
that  it  is  now  accepted  as  a  parameter  representing  the  nature  of 
the  surface  over  which  the  wind  (or  other  fluids)  is  flowing. 

In  fact,  tables  of  zq  (surface  roughness)  have  been  prepared  from 
application  (usually  erroneously)  of  the  adiabatic  profile  law, 
for  various  types  of  surfaces. 

Extension  of  eddy  coefficients  determination  above  the  layer 
of  constant  (with  elevation)  flux  values  is  also  being  pursued. 
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The  basic  relationship  is  a  modification  of  the  equation  suggested 
by  Clayton,  Covey,  and  Merryman  (1957) 


m,  z 


_ _  z  < .  Z\ 2 
C  \,l  7(1  “  Z> 


16) 


where  Z  is  the  vertical  extent  of  the  frictional  boundary  layer, 

l  is  the  vertical  extent  of  the  constant  flux  layer,  and  t  <  z  <  Z. 

Use  of  16)  has  been  limited  until  now  by  too  close  coupling 

of  K  with  K  p  and  unrealistic  values  of  K  «.  The  latter  has 
m,z  m,<C 

been  resolved  by  the  modified  layer  Richardson  number  and  improvement 
of  the  former  is  being  sought  through  modification  of  16)  by  a 
stability  parameter  as  indicated  below. 

K  „  "  c  Km  /  7(1  "  "  YRi)±n  17) 

In  this  equation,  n  is  an  integer  or  an  integer  reciprocal. 

Ri  is  a  modified  layer  Richardson  number  of  the  form  given  in  15) . 

Investigation  in  this  area  has  not  progressed  sufficiently  to 
justify  conclusions  as  to  useful  applicability. 

An  approach  to  K  determinations  in  the  frictional  boundary 

Itl 

layer  for  wider  extremes  of  stability  is  suggested  by  Professor 
Deacon.  This  approach,  in  part  empirical,  is  directly  addressed 
to  the  Low  Level  Meteorological  Model. 

Direct  measurements  of  shearing  stress  as  shown  in  Fig.  V-l  as 


well  as  data  from  other  researchers  (e.g.  Deacon  1953)  support  the 
following  relationship 


kf(Ri) 


18) 


z3S/3z  m 

Results  at  more  unstable  conditions  than  Ri  ■  -0.4  are  as  yet 
lacking  as  the  usual  unsteadiness  of  strongly  unstable  conditions 
makes  It  difficult  to  secure  reliable  observations.  On  the  stable 
side,  beyond  Ri  -  0.1,  turbulence  decreases  further  and  recent 
work  (Okamoto  and  Webb,  in  press)  shows  that  turbulence  in  the  usual 
sense  has  decreased  to  zero  at  a  height  of  2  m  when  Ri  ■  0.2. 
Radiative  heat  exchange  is  important  when  turbulence  is  small  so 
18)  is  not  expected  to  hold  beyond  about  Ri  -  0.15. 

Some  simple  linear  approximations  to  f(Ri),  as  noted  previously, 
have  been  used  for  near-neutral  conditions  but  if  applied  over  a 
wider  stability  range  are  found  to  be  unsatisfactory.  It  is  there¬ 
fore  feasible  to  investigate  a  non-linear  form  despite  the  added 
complexity  in  computation.  Further,  no  formulation  has  been 
proposed  which  is  valid  for  both  stable  and  unstable  conditions 
and  it  is  necessary  to  treat  them  separately. 

A  formula  proposed  by  Webb  for  unstable  conditions  fits  the 
data  of  Fig.  V-l  quite  well  but  the  'Keyps*  formula  is  but  little 
different  as  may  be  seen  from  Lumley  and  Panofsky  (1964,  p  113, 

Fig.  3.3)  and  the  latter  has  been  used  herein.  This  is 

^  -  k(l  -  Y'Ri)1"  19) 

in  which  y*  ■  18  is  suggested  by  Panofsky  on  the  basis  of 
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observational  material  at  O'Neill  and  elsewhere.  Equation  19)  can 
also  be  written  as 


p4  +  Y'(fr)P3 


20) 


in  which  P  - 

u* 

u*(3S/5z)6 
■"d  L'  ■  k8 

Oti  the  assumption  that  the  ratio  of  the  eddy  coefficients  of  heat 
and  momentum  is  constant  (y)»  then  L'  is  related  to  the  Obukhov 
L  by 

L'  -  yL  21) 

In  the  following,  use  is  made  of  the  integrated  form 

u. 

S  m  —  'J'(p-) ) 

o 

The  function  <i(z/T.’)  has  been  computed  by  Panofsky  and  is  reproduced 
in  Fig.  V-2. 

An  approximate  solution  for  20)  for  (z/L')  <  1  is  given  by  7) 
which  is  assumed  valid  for  stable  conditions  with  a  ■  5. 

With  other  than  linear  approximations  to  18)  there  is  the 
difficulty  that  explicit  relationships  for  and  K  are  unobtainable. 
This  arises  from  the  fact  that  u*,  for  a  given  Sg,  depends  on  the 
stability  conditions  from  the  surface  up  to  8  m  and  the  stability 
conditions  depend  not  only  on  temperature  gradient  but  also  on  the 


value  of  uA.  This  makes  it  necessary  to  proceed  by  aucceaaive 

approximations  as  will  be  seen  later. 

Utilising  the  preceding,  the  problem  is  to  find  u.,  K  «  and 

w  m,o 

Dg  given: 

Sg  -  wind  speed  at  8  m  height 

0g,  6^  ■  potential  temperatures  at  8  m  and  at  the  surface,  and 
sq  ■  the  roughness  parameter  of  the  surface. 

The  stability  criterion  to  be  employed  is  a  layer  Richardson 
number  for  the  surface  to  8  m  layer  as  shown  in  1A).  With 
similarity  between  wind  and  temperature  profiles  it  is  easily  shown 
that  8 

<‘/L'>8  '  ^  »> 

The  procedure  in  establishing  the  relationship  between  u^/Sg 
8 

and  R1  is  as  follows' 
o 

3 

1)  A  rough  estimate  of  u./Sfl  is  made  for  a  given  value  of  Ri 

w  o  o 

to  obtain  (z/L')g  from  Equation  23) 

2)  The  (z/L')g  estimate  gives  <J>g  from  Fig.  V-2 

3)  With  the  given  z  for  the  surface  and  the  value  of  u./S_ 

O  o  w  o 

is  obtained  from  Equation 

A)  The  improved  estimate  of  uA/Sg  is  employed  in  Steps  1  to  3 
until  no  significant  change  occurs 

For  zq  ■  2  cm  the  resulting  relationship  is  the  full  line  of  Fig.  V-3. 

To  obtain  K  knowing  u*  one  makes  use  of  the  flux  equation 
in 

u  2  -  K  3S/3z 
*  m 
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which  gives 


,V  ,  U*  ,  . 

^Sg^zas/sz^a 


24) 


A  graph  of  P  versus  z/L'  constructed  from  Equation  20)  gives  the 

value  k/P  which  is  the  quantity  in  the  second  bracket  of  Equation  24). 

A a  u./S_  has  already  been  obtained,  K  0/S0  is  found.  For  z  •  2  cm 
*  o  m,o  ©  o 

the  resulting  relationship  is  the  full  line  of  Pig.  V-4. 

The  exchange  coefficient  Dg  defined  by 

“.2  -  »«*«  25) 

is  obtained  from  u^/Sg  as 


<?>2 

b8 


26) 


This  graphed  in  Fig.  V-5. 

To  apply  these  results  in  the  LLMM  analog  computations  it  is 

necessary  to  fit  simple  mathematical  expressions  to  the  curves  of 

Fig.  V-3  and  V-4.  Expontential  expressions  appear  the  most 

suitable  if  different  expressions  are  used  for  positive  and  negative 
8 

Ri.  They  can  be  chosen  to  make  u*/Se  and  K  „/Sa  approach  reasonable 

O  “OTOfOO 

8 

asymptotic  values  at  large  negative  Ri,  An  exponential  decrease 

8 

of  these  values  with  increasing  positive  Ri  also  has  the  merit  of 
avoiding  the  negative  values  which  can  occur  using  linearized 
relationships . 

The  relationships  shown  in  Fig.  V-3,  V-4,  and  V-5  as  broken 
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lines  are  listed  below 


Unstable  conditions  and  2  ■  2  cm 

o 

8 

u./S.  -  0.0885  -  0.020  exp(13.5  Ri)  27) 

w  o  o 

D0/S0  -  {8.05  -  3.35  exp  (9.?  Ri)}  x  10-3  28) 

00  o 

8 

K  _/S0  -  92.2  -  70  exp <7. 15  Ri)  29) 

in,  o  o  o 

Stable  conditions  and  z  ■  2  m 

o 

8 

u*/S-  «  0.0685  exp (-5.0  Ri)  30) 

*  o  O 

8 

Dg/Sg  -  0.00470  exp (-10.0  Ri)  31) 

K  J S_  -  22.2  exp(-25.0  R?)  32) 

m,<5  c  o 

Similar  calculations  can  be  made  for  a  number  of  other 

roughness  values  and  the  coefficients  in  Equations  27)  to  32) 

graphed  against  2  . 

o 

Table  V-5  shows  the  various  ratios  to  the  8  m  wind  for  z  ■  2  cm. 


Table  V-5. 

Ratios  of  u.f  K  «. 

w  in,  0 

for  z  ■  2  cm 

0 

and  Dg  to 

the  8  m  Wind  Speed 

Computed 

8 

Ri 

0 

<*/L')g 

u* 

S8 

S8 

°8  _3 

—  x  10 

S8 

+.10 

1.202 

.0341 

1.6 

1.16 

+.06 

0.515 

.0478 

4.4 

2.28 

+.04 

0.300 

.0547 

7.2 

2.99 

+.02 

0.133 

.0618 

12.0 

3.82 

0 

0 

.0685 

22.5 

4.69 

-.01 

-0.058 

.0710 

28.5 

5.04 

CN 

O 

• 

1 

-0,112 

.0730 

33.5 

5.33 

1 

O 

-0.165 

.0745 

37.5 

5.55 
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Table  V-5 — Continued 


8 

Ri 

o 

(z/L‘ )g 

u* 

S8 

^8 

S8 

D8  .n3 
■=—  x  10 

S8 

-.05 

-  .266 

.0771 

44.0 

5.94 

-.07 

-  .362 

.0794 

50.0 

6.30 

-.10 

-  .500 

.0820 

57.0 

6.72 

-.15 

-  .715 

.0858 

67.0 

7.36 

-.20 

-  .910 

.0895 

75.0 

8.01 

Direct  comparison  of  Tables  V-4  and  V-5  cannot  be  made  since 

the  range  of  layer  Richardson  numbers  In  the  former  Is  limited. 

Approximately,  however,  K  „  values  from  Table  V-5  are  nearly  a 

m,  o 

factor  of  2  greater  for  negative  Ri  than  those  of  Table  V-4,  an 
order  of  magnitude  less  for  positive  Ri,  and,  of  course,  good 
agreement  exists  fo.-  Ri  ■  0. 

The  question  of  an  acceptable  diabatic  ind  profile,  for  all 
stabilities,  is  far  from  settled  and  investigations  continue. 

An  empirical  approach  based  on  data  collected  at  Stations  A  and  B 
during  their  Dallas  employment  is  the  subject  of  a  master's  thesis 
by  Mr.  Richard  Williams.  It  is  expected  that  this  thesis  investi¬ 
gation  will  be  completed  in  January  1969. 

Several  other  areas  are  under  study  in  the  analog  modeling 
activities.  Assessment  of  these  efforts  will  be  deferred  pending 
real-time  forecast  evaluations  of  the  LLMM  which  are  scheduled  to 
begin  In  early  1969.  These  efforts,  only  partially  under  Themis 
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sponsorship,  will  be  done  cooperatively  with  the  Meteorological 
section  of  Dugway  Proving  Ground. 
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GLOSSARY 


b 

c 


8 

h 

H 

k 

K 

D 

L 

P 

Ri 

S 


a 

6 

Y 


turbulent  anomaly  to  mean  wind 

parameter  in  the  equation  for 

K  as  a  function  of  elevation 
m 

specific  heat  of  air  at  constant  pressure 

integral  exchange  coefficient  for  momentum 
between  surface  and  8  m 

gravitational  acceleration 

elevation  of  constant  flux  level 

convective  heat  flux 

Von  Karman  constant  (.4) 

eddy  viscosity 

turbulent  scale  length 

normalized  logarithmic  wind  shear 

layer  Richardson  number 

mean  wind  speed 

frictional  wind  speed 

elevation 

surface  roughness  length 

elevation  of  upper  limit  of 
frictional  boundary  layer 

parameter  in  log  plus  linear  diabatic 
wind  profile 

parameter  in  Deacon  wind  profile 

parameter  in  the  equation  for  K  as 
function  of  elevation  and  stabiTity 


(era/sec) 

n  on-dimens lonal 
(cal/gm  deg) 

(cm/sec) 

2 

(cm/sec  ) 
cm 

2 

(cel/cm  sec) 
non-dimensional 
(cm2 /sec) 
cm 

non -dimensional 
non-d  irrensional 
(cm/sec) 

(cm/sec) 

cm 

cm 

cm 

non-dimensional 

non-dimensional 

non-dimensional 
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deg  C 


mean  potential  temperature 
friction  temperature 
mean  potential  temperature  of  layer 
air  density 

horizontal  shear  stress 
stability  function 


deg  C 
deg  K 
(gm/cra3) 

2 

(dynes /cm  ) 
non-dimension el 
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VI.  TURBULENT  FLOW  DEFINITION  THROUGH  REYNOLDS  STRESSES  MEASUREMENTS 

William  H.  Clayton 

The  design  of  the  mobile  Reynolds  Stress  facility  has  been 
completed  and  construction  of  the  station  has  begun.  Initially,  the 
planning  called  for  pulse  sonics,  using  three  four-head  crossed- 
path  systems  at  each  of  two  measuring  levels.  The  crosaed-path 
systems  were  subsequently  abandoned  in  favor  of  a  single  transmitter 
with  six  receivers  arranged  in  a  three  dimensional  orthogonal 
pattern  at  each  level.  Inasmuch  as  choice  of  the  pulse  type  was 
largely  based  on  construction  of  a  prototype,  it  was  felt  that  a 
more  thorough  study  cf  continuous  wave  types  should  be  pursued, 
since,  the  latter  could  lead  to  somewhat  Bimpler  analyzing  and 
processing  circuitry.  It  was  later  found,  more  or  less  as  expected, 
that  the  continuous  wave  sonic  continued  to  introduce  excessive 
design  complications  to  eliminate  reflection  interference,  and 
consideration  for  this  type  sonic  was  abandoned  and  final  choice 
was  made  of  a  pulse  sonic  with  the  single  transmitter  and  six 
receivers,  per  level,  configuration. 

Analyses  of  various  possible  control  systems,  that  is,  small 
digital  computers,  was  quickly  completed,  and  twelve  different 
manufacturers  of  such  equipment  were  initially  considered.  However, 
design  criteria,  costing,  power  consumption,  and  size  reduced  the 
computer  choice  to  four.  These  were: 
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1)  Sigma  2-Scientific  Data  Systems 

2)  DEP  516-Honeywell 

3)  SEL  810A-Systems  Engineering  Laboratories 

4)  2116B-Hewlett  Packard 

All  of  these  computers  are  so-called  third  generation  types, 
employing  Integrated  circuitry  and  any  of  them  would  be  adequate 
for  the  required  purposes.  Bid  requests  were  submitted,  but 
then  deferred  when  it  was  realised  that  the  interfacing  problem 
was  too  heavily  dependent  on  the  computer  to  be  employed  to  permit 
feasible  project  construction,  and  it  was  decided  to  incorporate 
external  counting  registers  and  computer  interfacing  as  a  part  of 
the  computer  requirements.  Modified  specifications  were  submitted 
to  the  companies  noted  above  and,  after  some  reluctance  on  the 
part  of  one  or  two  of  the  manufacturers,  formal  bid  responses  were 
received.  The  general  specifications  for  the  control  computer 
were  as  follows. 

1.  Word  length  -  16  bit 

2.  Maximum  memory  cycle  time  -  3u  sec. 

3.  Memory  capacity  -  8K  (16K  as  option  if  available 
in  basic  machine  housing) 

4.  Hardware  multiply  and  divide 

5.  ASR-33  teletype  input  -  output  with  necessary  interface 
(also  ASR-35  -v  option) 

6.  Power  requireter.ts  -  115/230  VAC,  60  cycle,  single  phase 

7.  Fortran  IV  -  wage  requirements  as  additional  cost, 
hardware  or  software 

8.  Special  features  (in-board  or  out-board) 

a.  One  megacycle  clock 

b.  Sixteen  registers  (16  bit)  which  can  be  separately 
interrogated  and  cleared  by  the  computer  and  which 
will  separately  accumulate  clock  pulses  until  a 
stop  pulse  is  received  from  an  external  transducer. 
Each  register  must  be  supplied  with  an  overflow 
indicator  which  can  be  recognized  by  the  computer 
to  identify  incorrect  data. 
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The  successful  bidder  was  Che  Honeywell  Company ,  and  Che  process 
controller,  complete  with  exCernal  registers  and  interfacing, 
with  ASR  35  teletype  Input/output  was  ordered  in  November  of  this 
year. 

The  use  of  sixteen  registers  is  to  accomodate  the  twelve 
sonic  signals  (six  at  two  levels  of  measurement) ,  two  temperature 
signals  at  each  level,  and  two  spares. 

Temperature  measurements  will  be  made  with  fast  response 
thermistors  whose  outputs  are  converted  to  frequencies.  The 
measurement  of  temperature,  independent  cf  the  sonic  system,  is 
considered  important  inasmuch  as  the  sonic  temperatures,  obtained 
from  the  sonic  system,  are  dependent  on  vapot  content  of  the  eir. 

In  addition,  independent  temperature  measurements  will  permit 
heat  flux  computations,  using  the  two  spare  channels  for  vertical 
velocity  determlnatlone  by  other  censors,  from  two  unrelated  sour  ms. 

Planned  housing  for  the  mobile  facility  la  shown  in  fig.  VI  -l. 
This  housing  is  a  renovated  and  modified  form  of  e  former  AM/CKC- 
radio  station.  Initial  planning  had  called  for  the  use  of  a  hous  Uvjt 
of  the  pickup  camper  type,  but  the  no-cost  availability  of  ths  rati; 
housing  plus  additional  space  made  this  s  feasible  substitution. 

Station  mobility  is  obtained  through  mounting  of  the  housln j 
on  a  1  1/2-ton  truck  (obtained  through  the  General  Services 
Administration) . 

Planning  is  not  complete  as  to  tower  and  Instrument  support 
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construction  but  tentative  planning  calls  for  a  trailer  configuration 
similar  to  the  Rohn  HD1-3-54G.  This,  or  a  similar  type  of  tower, 
will  permit  reasonable  mobility  with  regards  to  station  transfer 
as  well  as  rapid  change  of  measurement  levels. 

The  sonic  heads  and  temperature  sensors  to  be  employed  in 
this  facility  are  still  under  study  and  final  selection  will  not 
be  made  without  consultation  with  other  investigators  in  the  sonic 
field.  In  this  regard,  visits  are  planned  in  early  1969,  to  the 
University  of  Washington  (Dr.  J.  A.  Businger),  Iowa  State  University 
(Dr.  Robert  Stewart  and  Dr.  Robert  Post)  and  to  the  Air  Force 
Cambridge  Research  Laboratories,  Hanscom  Field,  Massachusetts 
(Dr.  Chandran  Kainal) . 

The  preceding  discussion  plus  the  'ahead-of-schedule'  status 
of  this  Sub-task  is  not  intended  to  imply  that  completion  of 
this  facility  is  largely  mechanistic  and  straightforward.  Several 
problems  remain  to  be  solved  of  an  anticipated  nature  and  un¬ 
doubtedly  several  unanticipated  ones  will  be  met.  There  can  be 
no  question  that  a  large  part  of  the  difficulties  will  be  engine¬ 
ering  in  nature,  however,  the  theoretical  implications  of  sonic 
anemometry  are  at  present  only  unproved  postulates  and  evaluation 
of  these  postulates  is  not  expected  to  be  routine.  This  is  perhaps 
best  indicated  by  noting  that  all  of  the  researchers  noted  above 
have  constructed  successful  prototype  sonic  anemometers,  yet  none 
of  these  devices  have  been  consistently  and  gainfully  employed 
in  the  field. 
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VII.  MISCELLANEOUS 


In  the  preceding  review  of  Themis  activities  during  the  first 
year,  attention  has  been  focused  on  results  and  avoidance  of 
detailed  descriptions  of  unsuccessful  and/or  essentially  logistical 
endeavors.  This  is  particularly  true  for  the  data  collection 
activities  under  Sub-task  b,  the  calibration  and  installation  proce¬ 
dures  of  the  automatic  meteorological  stations  under  Sub-task  c, 
and  routine  model  evaluation  and  efforts  also  under  Sub-task  c. 

Some  review  of  these  activities  have  been  given  in  the  Quarterly 
Progress  Reports  and  omission  herein  is  strictly  based  on  editorial 
concept.  An  idea  of  the  magnitude  of  these  omitted  efforts  may  be 
afforded  by  noting  that  all  of  the  originally  progosed  efforts  are 
closely  following  the  proposed  schedule  with  the  deviation,  for 
the  most  part,  being  one  of  acceleration  rather  than  retardation. 

Two  new  efforts  will  be  initiated  during  the  second  year.  One 
under  Sub- task  b,  the  other  under  a  new  Sub- task  e. 

The  first  of  these  efforts  "Anomalous  Wintertime  Precipitation 
in  the  Western  United  States  and  Related  Meteorological  Variables", 
will  be  directed  by  Dr.  R.  A.  Clark.  The  research  itself  will 
largely  be  conducted  by  Mr.  L.  G.  Cobb  as  a  Ph.D.  dissertation.  A 
brief  resume  of  the  procedures  that  will  be  followed  in  conducting 
this  research,  given  below,  has  been  extracted  from  a  proposal 
submitted  by  Mr.  Cobb  to  the  Technical  Director. 
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The  proposed  research  will  Investigate  the  deviations  from 
the  climatic  mean  monthly  precipitation  of  the  wet  season  for 
the  Pacific  coast  of  the  United  States.  Relationships  between 
the  rainfall  anomalies  in  the  total  problem  will  be  investigated. 

The  monthly  Climatological  Data  National  Summary  contains 
maps  of  the  total  precipitation,  percentage  of  normal  precipi¬ 
tation,  average  sea-level  pressure,  and  departure  of  average 
pressure  from  normal,  along  with  the  average  height  and  tempera¬ 
ture  and  resultant  winds  for  standard  lsobaric  surfaces.  These 
will  be  used  to  determine  the  months  of  significant  anomalies 
(less  chan  75  percent  or  greater  than  150  percent  of  the  mean) 
of  rainfall  In  the  states  of  California,  Oregon,  and  Washington 
for  the  period  1952-1967.  The  wet  season  (October-April)  will 
be  considered  because  only  during  those  months  does  enough 
precipitation  occur  so  that  deviations  from  the  average  have 
real  meaning. 

From  monthly  mean  radiosonde  data  for  12  stations  located 
in  tha  western  United  States,  the  flux  of  water  vapor  across 
the  Pacific  coast  will  be  calculated  for  40  months  selected 
from  the  period  1952-1967.  Ten  months  which  are  representative 
of  each  -A  ‘he  following  four  situations  will  be  chosen  for 
analysis:  above-normal  precipitation  over  the  whole  west 
coast  region,  below-normal  precipitation  over  the  whole  west 
coast  region,  much-below-normal  precipitation  in  the  southern 
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half  of  the  region  and  normal  or  greater  precipitation  in  the 
northern  half,  and  much-below-normal  precipitation  in  the 
northern  half  and  above-normal  in  the  southern  half  of  the 
area  being  considered.  Calculations  of  the  flux  of  vecer 
vapor  will  be  made  for  the  surface  and  each  50-mb  level  above 
to  400  mb  by  the  technique  used  by  Cummings  and  Franceschlnl. 

It  is  expected  that  a  careful  analysis  of  this  data  will  reveal 
characteristic  features  of  the  flux  of  moisture  for  each 
station  and  situation,  and  will  shed  light  on  the  question  of 
whether  increased  precipitation  is  caused  by  a  more  efficient 
precipitation  mechanism  or  a  more  abundant  moisture  supply 
(or  both). 

Mean  monthly  Bea-surface  temperatures  are  available  for 
the  eastern  North  Pacific  Ocean  for  the  period  being  investi¬ 
gated.  These  will  be  studied  in  conjunction  with  the  analyses 
of  water-vapor  flux  to  see  if  the  water  temperatures  significantly 
affect  the  water  content  of  the  air  and  the  atmospheric  pressure 
patterns. 

The  rainfall  anomalies  will  be  statistically  correlated 
with  the  anomalies  of  the  average  sea-level  pressure  of  that 
region  to  det^hnine  conditions  associated  with  the  pra^  re 
field  which  are  conducive  to  abnormal  precipitation.  Most 
previous  investigators  have  attempted  to  relate  precipitation 
and  upper-air  patterns  (mostly  700-mb);  however,  the  present 


study  will  utilise  sea-level  pressure  because  the  pressure 
patterns  at  the  surface  reflect  the  behaviur  of  the  total 
atmosphere.  A  preliminary  Investigation  has  revealed  very 
good  correspondence  between  the  anomalous  component  of  the 
wind  flow  and  areas  of  abnormal  precipitation.  Once  the  patterns 
of  pressure  anomaly  that  are  associated  with  abnormal  rainfall 
have  been  identified,  theoretical  relationships  will  be 
utilised  to  show  what  is  occurring  in  the  overlying  atmosphere 
in  terms  of  changes  in  moisture  advection,  convergence,  and 
relative  vorticity.  The  results  of  this  study  will  be  compared 
to  previous  studies  of  a  more  general  nature. 

The  second  new  effort,  directed  by  A.  H.  Thompson,  is  titled 
"Tne  Change  of  Intensity  in  Cloud  Systems  Associated  with  Easterly 
Perturbations  as  They  Move  Over  Areas  of  Differing  Water  Temperature". 
The  research  will  be  largely  conducted  by  Mr.  J.  E.  Arnold  as  a 
Ph.D,  dissertation  subject.  A  brief  resume  of  the  procedures  to  be 
employed  in  conducting  this  research,  as  given  below,  have  been 
extracted  from  a  proposal  submitted  by  Mr.  Arnold  to  the  Technical 
Director. 

The  temperature  field  of  the  Caribbean  Sea  and  equatorial 
western  Atlantic  will  be  examined  through  the  use  of  commercial 
ship  data  available  from  the  National  Weather  Records  Center, 
Asheville,  N.  C. ,  and  from  conventionally  prepared  northern 
hemisphere  weather  charts.  It  has  been  found  that  commerical 


135 


ship  data,  when  averaged  over  a  7-day  period,  provide  a 
reasonably  good  representation  of  ocean  surface  temperature. 

From  composite  sea-surface  temperature  charts  of  this  nature  the 
temperature  of  the  ocean  surface  can  be  acquired  for  periods 
when  active  systems  are  propagating  through  the  area  of  interest 
It  should  he  pointed  out  that  although  the  temperature  of  the 
ocean  acquired  by  this  method  is  not  always  correct,  i.e.,  it 
may  be  too  warm  or  too  cold,  the  temperature  patterns  are 
comparable  to  those  obtained  by  research  vessels. 

The  vertical  thermal  structure  of  the  upper  layer  of  the 
ocean  also  will  be  considered  by  the  use  of  bathythermographs 
available  from  time  to  time  in  the  area  of  interest  and  from 
historical  data  on  the  vertical  thermal  structure  of  the  ocean. 
This  information  will  be  used  to  evaluate  the  capability  of 
the  ocean  to  supply  or  absorb  heat  during  the  periods  studied 
and  the  effects  that  wind  mixing  would  have  on  the  surface 
layers  on  the  oceans. 

Surveillance  of  cloud  systems  moving  through  the  area  of 
interest  will  be  accomplished  using  ESSA  and  Nimbus  cloud 
photographs.  Individual  disturbances  will  be  followed  and 
the  cloud  activity  recorded  in  terms  of  area  ard  type  of 
activity.  The  distribution  of  cloud  activity  about  the 
disturbance  will  be  examined  on  an  interdiurnal  basis  from  ESSA 
and  Nimbus  photographs  and  the  cloud  range  during  short  periods 
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i.e. ,  during  the  day,  will  be  studied  from  the  ATS  III  photo¬ 
graphs.  This  should  provide  informative  data  as  to  cloud 
distribution  about  a  system  over  a  water  area  of  uniform 
temperature  and  the  response  of  the  cloud  distribution  as  the 
disturbance  passes  over  a  region  of  strong  surface  temperature 
gradient.  In  addition  to  cloud  examination,  the  ATS  III 
photographs  will  be  used  to  determine  cloud  motion  in  the 
vicinity  of  the  disturbances  studied. 

A  composite  picture  of  the  atmosphere  will  be  obtained 
from  existing  radiosonde  and  surface  stations  situated  around 
the  region  to  be  investigated.  The  patterns  of  divergence 
obtained  from  the  rawindsonde  and  plbal  measurements  surrounding 
the  Caribbean  will  be  integrated  with  the  patterns  of  divergence 
obtained  from  cloud  motion  (assumed  to  be  representative  of  the 
wind  field  at  cloud  level)  to  arrive  at  an  estimated  field  of 
vertical  motion  in  the  vicinity  of  the  disturbance.  The 
derived  field  of  motion  in  the  atmosphere  and  the  sea-surface 
temperature  field  will  be  used  as  input  data  from  which  the 
life  of  the  disturbance  and  associated  cloud  intensity  will 
be  evaluated. 

With  the  above  data  it  is  proposed  to  take  a  period  when 
traveling  perturbations  in  the  easterlies  are  known  to  exist  in 
the  Caribbean  and  western  Atlantic  region  and  maintain  close 
surveillance  of  the  region  to  find  weak  or  moderate  disturbances 
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•ntarlng  Che  area  of  Interest.  It  is  anticipated  that  a 
sufficient  number  of  disturbances  will  be  found  in  one  season, 
June-October,  1968,  to  enable  an  evaluation  of  the  proposed 
technique  to  determine  the  relationship  between  cloud  activity. 
Intensity  of  the  system,  and  sea-surface  temperature  field 
over  which  the  perturbation  passes. 

The  western  equatorial  Atlantic  and  the  Caribbean  Sea 
have  been  chosen  as  the  geographic  region  for  Investigation 
because  of  the  availability  of  data.  The  Caribbean  Sea  Is 
surrounded  by  rawlnsonde  and  surface  reporting  stations, 
enabling  a  basic  framework  of  upper-air  and  surface  stations 
to  be  used  in  connection  with  the  satellite  data.  In  addition, 
there  is  sufficient  shipping  within  the  area  to  enable  reasonably 
accurate  evaluation  of  the  sea-surface  temperature  field.  The 
region  also  la  near  the  subpoint  of  the  ATS  III  satellite,  a 
feature  definitely  advantageous  In  determining  cloud  motion 
from  the  (Image  Dissector  Camera  System)  pictures.  It  is 
generally  accepted  also  that  the  area  of  the  Caribbean  Sea 
is  within  the  normal  tracks  of  perturbations  as  the  systems 
pass  from  the  Atlantic  toward  the  Yucatan  Channel.  Preliminary 
investigation  of  the  Caribbean  sea-surface  temperatures  as 
determined  from  commercial  ship  dote  and  inspection  of  the 
oceanographic  atlases  reveal  that  sufficient  gradients  of 


138 


sea-surface  temperature  exist  through  the  proposed  time-period 
of  study  to  evaluate  realistically  the  proposed  hypothesis. 

It  is  intended  to  obtain  ATS  III  photographs  from  the 
National  Aeronautics  and  Space  Administration  for  the  construction 
of  film  strip  loops  to  determine  cloud  motion.  Tentative 
arrangements  have  been  made  with  Mr.  L.  F.  Hubert,  ESSA,  to  have 
the  National  Environmental  Satellite  Center  construct  the  film 
strip  loops  and  send  copies  to  Texas  A&M.  A  loop  projector 
system  similar  to  that  developed  by  Fujita  at  the  University 
of  Chicago  will  be  constructed  at  Texas  A&M.  The  technique 
of  analysis  described  by  Fujita  will  be  used  tc  determine  . 
cloud  motion. 

Finally,  it  should  be  noted  that  the  activities  of  the  second 
year  will  follow  proposed  schedules  but  with  somewhat  less  rigidity 
than  that  employed  during  the  first  year.  The  efforts  of  the  first 
year  have,  In  a  sense,  been  exploratory  and,  as  anticipated,  new 
lines  of  research  have  been  opened,  and  these  will  be  followed  as 
time,  staff,  and  facilities  will  allow.  Primarily,  these  are 
changes  in  the  originally  proposed  procedures  which,  necessarily, 
were  tentative,  and  no  basic  changes  in  direction  of  the  Themis 
efforts  at  Texas  A&M  are  currently  planned. 
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